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An  improved  cross-shore  beach  profile  evolution  model  was  developed  based  on 
equilibrium  beach  profile  and  surface  roller  concepts.  A more  physically  sound  wave 
transformation  model  and  roller  model  were  adopted  to  improve  the  calculation  of  wave 
height  transformation  and  energy  dissipation  inside  the  surf  zone.  Including  the  surface 
roller  model  produced  a new  equilibrium  beach  profile  shape  which  can  represent  bar 
profiles  to  a limited  degree.  The  profile  scale  factor  in  the  equilibrium  beach  profile 
formulae  was  modified  through  the  fall  velocity  parameter  based  on  best  fits  of  measured 
profiles  from  large-scale  laboratory  experiments.  To  predict  beach  profile  response  under 
wave  action,  a finite  difference  method  was  applied  to  solve  the  sediment  transport  and 
continuity  equations  numerically  through  an  implicit,  double  sweep  scheme.  After  model 
calibration  with  eleven  large-scale  wave  flume  experimental  data,  it  was  found  that  a 
nonlinear  transport  relationship  with  n = 2 provided  better  overall  predictions  than  n 
values  of  1 and  3.  This  nonlinear  relationship  provided  a reasonable  explanation  for 


different  lime  scales  of  beach  profile  evolution  under  different  wave  and  initial  profile 
conditions. 

The  model  was  applied  to  additional  large-scale  wave  flume  experimental  data 
from  SUPERTANK  and  LIP  1 1 D and  field  data  afler  major  storms  ut  Ocean  City, 
Maryland  and  the  Panhandle  Coast  of  Florida,  along  with  comparisons  with  two  other 
closed  loop  models  (SBEACH  and  CROSS).  Results  showed  that  all  three  models  had 
different  degrees  of  underprediction  of  volume  changes  and  contour  recessions  afler 
storms.  The  model  can  generate  bar  profiles  under  both  regular  and  random  wave 
conditions,  which  is  similar  to  the  SBEACH  model  but  no  empirical  assumption  is  made 
in  the  wave  transition  region. 


CHAPTER  I 

INTRODUCTION  AND  OBJECTIVE 
I.I  Introduction 

Sand  beaches  are  valuable  resources.  They  not  only  provide  protection  front  severe 
storms  and  hurricanes,  but  also  generate  substantial  revenue  as  a recreational  amenity. 
Under  the  continuous  seasonal  action  of  waves,  currents  and  winds,  beaches  often  erode 
during  winter  and  accrete  during  summer.  The  sediment  transport  at  a point  in  the 
nearshore  zone  includes  cross-shore  and  longshore  components.  Cross-shore  sediment 
transport  includes  both  offshore  transport  such  as  occurs  during  storms  and  onshore 
transport  during  mild  wave  climates.  Transport  in  these  two  directions  occurs  in 
significantly  distinct  modes  (suspended  load  dominates  offshore  transport,  while  bedload 
dominates  onshore  transport)  and  with  remarkably  disparate  time  scales.  The  offshore 
transport  moves  sand  from  the  upper  profile  to  offshore  in  hours  to  days  during  a storm, 
whereas  the  onshore  transport  may  take  months  to  years  to  move  sand  to  the  dry  beach  in 
which  a "ridge  and  runnel"  system  often  develops.  Difficulties  in  predicting  transport  in 
the  two  directions  also  differ  substantially.  Offshore  transport  is  the  more  predictable  of 
the  two.  This  is  fortunate  for  coastal  engineers,  since  there  is  considerably  greater 
engineering  concern  and  interest  in  offshore  transport  because  of  its  potential  damage  to 
structures  and  loss  of  land. 

Cross-shore  sediment  transport  is  relevant  to  a number  of  coastal  engineering 
problems  including  the  following: 


Cross-shore  equilibration  of  a beach  nourishment  project. 


• Profile  nourishment  in  which  the  sand  is  placed  in  the  nearshore  with  the 
expectation  that  it  will  move  landward  nourishing  the  beach. 

• Shoreline  response  to  sea  level  rise. 

• Seasonal  shoreline  changes. 

• Overwash,  the  process  of  landward  transport  due  to  overtopping  of  the  normal  land 
mass  due  to  high  tides  and  waves. 

• Scour  immediately  seaward  of  shore  parallel  structures,  such  as  seawalls. 

• Cross-shore  and  longshore  coupled  transport  around  shore-normal  structures  in 
which  the  steeper  and  milder  slopes  on  the  updrift  and  downdrift  sides  induce 
seaward  and  landward  transport  components,  respectively. 

One  successful  and  commonly  adopted  method  to  protect  storm  damage  and  attract 
more  people  to  the  beach  is  beach  nourishment,  especially  in  Florida.  Beach  nourishment 
projects  restore  the  eroded  beach  with  additions  of  large  quantities  of  good  quality  sand 
to  the  natural  system.  Beach  nourishment  profiles  arc  usually  constructed  with  slopes 
steeper  than  equilibrium.  The  time  scale  for  profile  equilibration  is  very  significant  to 
design  and  interpretation  of  beach  nourishment  project  revolution,  because  wider  beaches 
during  equilibration  provide  more  area  for  recreation  and  a greater  degree  of  storm 
protection  to  upland  structures.  Additionally  profile  equilibration  can  be  misinterpreted 
as  a loss  of  material  and  consequently  a poor  and  unanticipated  performance, 

A clear  understanding  of  the  time  scales  of  profile  equilibration  provides  guidance 
in  the  beach  nourishment  project  design  process.  Examples  include  the  volume  of  sand 
to  be  placed  in  order  to  attain  the  desired  beach  width  during  the  life  span  of  the  project, 
and  the  best  design  slope  of  the  placement  to  ensure  that  less  sand  is  carried  offshore. 

This  understanding  is  also  helpful  for  analysis  of  beach  nourishment  project  performance 
and  for  regulatory  purposes  such  as  establishing  coastal  construction  control  lines. 


Forces  governing  beach  profile  equilibration  processes  include  natural  and  artificial 


forces.  The  natural  forces  include  tides,  waves,  and  sea  level  rise.  The  most  significant 
artificial  force  is  due  to  profiles  placed  at  slopes  sleeper  than  equilibrium  that  induces 
seaward  sediment  transport  and  associated  shoreline  rctrcat.  Due  to  the  complexities  of 
surf  zone  hydrodynamics,  sediment  characteristics,  and  wave  conditions,  an  analytical 
treatment  of  beach  profile  change  is  very  difficult.  Thus,  numerical  modeling  (which 
allows  flexibility  in  the  specification  of  initial  and  boundary  conditions  and  the 
representation  of  arbitrary  forcing)  is  a powerful  and  appropriate  tool. 

1.2  Characteristics  of  Bench  Systems 
1.2.1  Approach  to  Equilibrium 

Beach  profiles  can  attain  equilibrium  states  under  constant  environment  conditions. 
An  equilibrium  profile,  by  definition,  is  approached  when  the  cross-shore  transport  rates 
at  all  points  become  zero,  as  arc  the  longshore  gradients  in  longshore  transport.  So  the 
decay  of  cross-shore  transport  rate  with  time  is  one  direct  measure  of  profile 
equilibration.  Figures  1-1  and  I -2  show  two  cases  of  transport  rate  distribution  with  time 
calculated  fiom  profile  measurements,  one  for  an  crosional  case  and  one  for  an 
accretionary  case  from  LVVT  experiments  (Kraus  and  Larson,  19SS).  Other  quantities 
that  can  be  used  to  represent  the  profile  equilibration  process  include  eroded  volume  (or 
bar  volume)  for  crosional  profiles  and  berm  volume  for  accretionary  profiles,  shoreline 
position  changes,  bar  positions,  and  so  on.  These  confirm  that  the  equilibrium  concept  is 
valid  and  that  a numerical  model  developed  for  simulating  realistic  beach  profile  change 
needs  to  include  this  property. 


1.2.2  Equilibrium  Beach  Profiles 

A quantitative  understanding  of  the  characteristics  of  equilibrium  beach  profiles  is 
central  to  rational  design  of  many  coastal  engineering  projects  and  to  the  interpretation  of 
nearshore  processes.  Several  features  of  equilibrium  profiles  are  well  known: 

• They  tend  to  be  concave  upward. 

• Smaller  and  larger  sand  diameters  are  associated  with  milder  and  steeper  slopes, 
respectively. 


The  beach  face  is  approximately  planar. 


• Sleep  waves  and  finer  sands  result  in  milder  slopes  and  a tendency  for  bar 

formation. 

1.2.2. 1 Statistical  approach  to  equilibrium  beach  profiles 

Through  analysis  of  beach  profiles  from  the  Danish  North  Sea  coast  and  Mission 
Bay.  California.  Bruun  (1954)  found  that  profiles  tended  to  follow  the  simple  relationship 

h = Ay’“  (1.1) 

where  A is  a dimensional  scale  parameter  with  units  of  length1 3 and  depends  primarily 
on  sediment  characteristics. 

Using  502  beach  profiles  from  the  Atlantic  and  Gulf  Mexico  shorelines  and  a least- 
square  procedure,  Dcan(1977)  fitted  an  equation  ofthe  form  h = Ay".  A central  value 
of  m = 2/3  was  found,  which,  according  to  linear  wave  theory,  is  consistent  with  uniform 
wave  energy  dissipation  per  unit  volume  within  the  surf  none  as  the  dominant  destructive 

Bodge  ( 1 992)  proposed  an  exponential  expression  for  equilibrium  profile  as 

A-«l-«-*)  (1.2) 

where  B and  k are  constant  parameters.  This  expression  along  with  the  2/3  power 
equilibrium  profile  were  fit  to  a set  of  data  from  504  beach  profiles  from  the  U.S.  east 
coast  and  Gulfof  Mexico  which  Dean  (1977)  originally  used  to  investigate  the  h = Ay1" 
shape.  Bodge  found  that  most  (80%  to  86%)  were  better  fit  by  the  exponential  form 
relative  to  the  2/3-powcr  shape.  He  concluded  that  the  improved  fit  of  the  exponential 
expression  was  likely  to  be  most  notable  for  beaches  with  decreasing  grain  size  in  the 
offshore  direction.  The  leading  coefficient  B may  yield  insight  to  the  depth  of  closure; 
and  the  coefficient  k may  describe  sediment  type  or  its  gradation  in  the  offshore 


direction.  However,  no  attempt  was  made  to  correlate  the  parameters  with  beach  or  wave 
characteristics. 

Komar  and  McDougal  ( 1 994)  considered  an  alternative  form  of  the  exponential 
beach  profiles  (Eq.1,3). 


where  So  is  the  beach  face  slope.  They  noted  that  because  there  have  been  many  past 
studies  that  relate  So  to  beach  sediment  grain  sizes  and  wave  conditions,  this  exponential 
beach  profile  form  only  contains  one  adjustable  coefficient,  k,  whose  evaluation  can  be 
based  on  the  coordinates  of  the  closure  depth  of  profile  change  or  equally  on  some 
arbitrary  chosen  offshore  depth  coordinate. 

By  reasoning  that  the  forcing  mechanism  landward  and  seaward  of  the  breakpoint 
bar  (by  shoaling  waves  and  breaking  waves  respectively)  can  be  significantly  different, 
Inman  ct  al.  (1993)  divided  the  beach  profile  into  two  independent  portions  separated  by 
the  bar  crest.  The  two  portions  were  termed  as  bar-berm,  from  the  berm  crest  to  the 
breakpoint  bar;  and  shorense.  seaward  of  the  breakpoint  bar  (Figure  1-3). 


0-3) 


Figure  1-3.  Definition  sketch  of  bar-berm  and  shorerise  profile  (Inman  el  al.,  1993). 


Two  power  functions  similar  to  h = Ay”  were  used  independently  to  reconstruct  the 
two  segments.  It  was  demonstrated  that  the  barred  profiles  were  reproduced  better, 
especially  in  the  vicinity  of  the  bar  and  the  offshore  portion,  than  the  one  segment  model 
by  h - Ay”.  This  profile  model  has  8 parameters  requiring  quantification. 

Wang  and  Davis(1998)  developed  a three-segment  beach  profile  model  based  on 
the  Inman  etal.  (1993)  approach  for  reproducing  the  commonly  observed  bar  and  trough 
features.  The  model  was  developed  under  the  general  form  of 

A=FW)  (1.4) 

where  the  function  F is  determined  empirically  based  on  regional  boundary  condition. 
The  three  segments  used  are 

• Inner  surf  zone,  from  the  lowest  elevation  of  the  trough  to  the  shoreline,  dominated 
by  secondary  wave  shoaling  and  breaking. 

h = A,y-,  0<y<y„  (1.5) 

• Landward  of  the  bar  crest,  characterized  by  a landward  dip  and  a landward  increasing 
depth.  The  dominant  hydrodynamic  process  is  rapid  energy  dissipation  at  the  bar  lop 
due  to  breaking  and  wave  reformation  caused  by  increasing  water  depth. 

h=li„*  (y-y„),  yu<y<ybi  (1.6) 

•V*  -y. 

• Nearshore,  from  the  bar  crest  to  the  closure  depth,  this  portion  is  dominated  by  wave 
shoaling  and  breaking  on  the  bar  lop. 

*-40--*r.  y«<y<y^  0-7) 

where  Ai  and  A;  are  dimensional  scale  parameters  related  to  grain  size  for  inner  surf  and 
nearshore,  respectively;  ini  and  m2  are  empirical  shape  parameters  controlling  the  beach 
slopes;  ha  and  yu  are  water  depth  at  trough  and  its  distance  to  the  shoreline;  hbi  and  ysi  arc 
water  depth  at  bar  top  and  its  distance  to  the  shoreline;  y„.  is  the  distance  from  shoreline 


to  the  closure  depth  of  the  profile  and  y<  is  the  intercept  of  the  nearshore  portion  with 


zero  water  level.  The  quantity  m,  was  found  to  be  well  approximated  by  2/3,  while  mi 
was  larger  than  2/3,  with  an  averaged  best-fit  value  of  0.82.  Application  of  this  model 
requires  knowledge  of  10  parameters,  none  of  which  can  be  determined  without 
availability  of  the  profile.  Thus,  this  model  can  be  applied  in  a diagnostic,  but  not  a 

1.2.2.2  Based  on  equilibrium  mechanisms 

Bowen  (1980)  established  an  equilibrium  beach  profile  model  using  an  energetic 
model  of  suspended  sediment  transport  based  on  wave  orbital  motions  and  a perturbation 
or  drift  velocity  within  the  nearshore.  The  first-order  solution,  assuming  symmetrical 
wave-orbital  motions,  yielded  a beach  profile  having  a 2/3-power  depth  variation. 
However,  when  the  higher  harmonic  terms  were  included  to  simulate  the  more  realistic 
asymmetrical  wave-orbital  motions  as  found  in  Stokes-type  waves,  the  analysis  then 
yielded  a profile  having  a 2/5-power  depth  variation. 

Based  on  Bagnold's  energetic  approach  for  bed  load  and  suspended  load  under 
steady  flows,  Leont'ev  (1985)  developed  an  equilibrium  profile  model  based  on  the 
following  considerations:  (I)  Bottom  sediment  grains  are  homogeneous  in  density  and 
size;  (2)  Bedload  sediment  transport  is  onshore  due  to  the  wave  orbital  velocity 
asymmetry;  (3)  Suspended-load  transport  is  offshore  and  its  rate  is  equal  to  the  product  of 
the  suspended  sediment  concentration  and  the  undertow  velocity.  The  model  results  in 
Eq.1.8. 


(1.8) 


where  y is  positive  onshore  with  origin  at  the  breaking  point,  h is  water  depth,  ho  is  the 
water  depth  at  wave-breaking  point,  y = H/h,  H is  wave  height,  w is  the  particle  fall 
velocity,  Co  is  the  wave  phase  speed  at  the  break  point  and  B is  a constant  of  the  older  of 
10"1,  Using  two  different  y forms  at  the  offshore  (Green's  law)  and  onshore  region 
(constant)  from  the  breaking  point,  an  equilibrium  profile  was  obtained  which  was 
convex  in  the  offshore  region  (where  wave  energy  grows)  and  concave  in  the  surf  zone 
(where  wave  energy  dissipation  occurs).  This  model  was  applied  to  both  monochromatic 
and  random  waves.  Figure  1 -4  compares  measured  and  theoretical  results. 


Figure  1-4,  Comparison  between  measured  and  theoretical  EBP  (Leont’ev,  1985).  Solid 
line  i and  dashed  line  2 refer  to  measured  and  calculated  profile,  respectively. 

Larson  (1988)  modified  the  concept  of  Dean’s  energy  dissipation  per  unit  volume 
by  replacing  the  spilling  breaker  assumption  with  the  breaking  model  by  Dally  et  al, 

( 1 984),  which  considers  the  dissipation  due  to  breaking  in  the  surf  zone  to  be  given  by 

dy  h 

where  the  energy  flux  F = ECt , F,  is  the  stable  energy  tlux  based  on  a stable  breaking 
wave  height  H , = IV) . Equating  this  energy  flux  relationship  to  Dean's  uniform  energy 
dissipation  per  unit  volume  leads  to  the  resulting  equilibrium  profile 


(1.10) 


r-ttr 


(1.11) 
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condilions  of  equilibrium  and  the  governing  mechanisms.  The  resulting  EBP  shapes 
under  breaking  and  nonbreaking  waves  are  given  by 

h = Aym,  0<y<yi  (1.14) 

A =[V" + *(>' -*))'.  y>y»  (115) 

where  A and  B are  dimensional  shape  parameters.  A is  related  to  grain  size  (or  fall 
speed);  whereas  B is  a function  of  offshore  wave  conditions  and  sediment  characteristics. 
The  values  of  the  powers  were  determined  to  be  m = 2/3  and  n = 0. 1 5 to  0.30.  Note  that 
the  theoretical  EBP  models  were  essentially  based  on  monochromatic  (or  representative) 
wave  condilions  and  the  effects  of  wave  randt  n s ere  not  addressed  explicitly. 

Medina  ct  al.  (2000)  proposed  a new  equilibrium  beach  profile  formulation  that 
treats  the  shoaling  portion  of  the  profile  independently  from  that  of  the  breaking  portion. 
The  formulation  considers  wave  reflection  from  the  beach  and  the  variation  of  sea  level 
due  to  tide.  The  theoretical  basis  of  the  model  is  the  time-averaged  wave  energy  flux 
equation  for  straight  and  parallel  bottom  contours. 


where  F is  the  net  shoreward  energy  flux  per  unit  width  and  e is  the  energy  dissipation 
rate  per  unit  surface  area.  The  equation  involves  three  variables:  wave  height,  water 
depth  and  wave  energy  dissipation.  Two  different  wave  energy  dissipation  models  and 
wave  height  variations  across  the  profile  are  proposed  for  breaking  and  non-breaking 
waves,  respectively. 

• EBPs  under  breaking  waves 


(1.17) 


H=rh 


<1.18) 


J,=(^)1"+^AI.°cy<>i. 


(1.19) 


(1.20) 


where  Db’  is  the  wave  energy  dissipation  per  unit  water  depth  caused  by  wave  breaking,  y 


is  the  brcakcr-to-dcpth  ratio,  B is  a constant  which  depends  on  wave  period,  and  is  related 
to  the  reflected  energy  flux  and  is  zero  if  reflection  is  negligible. 

• EBPs  under  nonbreaking  waves 


®=d;  (i.2i) 

H = pjh  (1.22) 


' y - p.  - fc)1" + . y^y'y. 


(1.23) 


(1.24) 


3'.=J,.0-^T)  + ^r*.‘  (1-25) 

where  Df  is  the  constant  wave  energy  dissipation  per  unit  area  caused  by  friction,  p and 
D are  constants.  Similar  to  B in  the  breaking  zone,  D is  also  related  to  reflected  energy 
flux  and  is  zero  if  no  reflection.  The  subscript  “a”  refers  to  the  offshore  limit  of  the 
shoaling  profile. 

An  increase  in  wave  height  will  have  two  effects  on  the  equilibrium  shape.  First, 
the  shoaling  profile  will  become  steeper  as  the  shape  parameter  C increases.  Second,  the 
offshore  limit  of  the  surf  profile  will  also  increase  due  to  the  greater  breaking  depth. 
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Besides  these,  wave  reflection  adds  a new  link  between  both  portions  of  the  profile.  A 
variation  in  the  wave  period  will  change  the  shape  parameters  B and  D and  consequently, 
the  overall  shape  of  the  profile.  Best-fit  tests  using  the  above  EBP  shape  against  over  50 
profi  les  from  1 3 beaches  along  the  Spanish  coast  were  conducted  to  detennine 
relationships  among  the  shape  parameters  (A.  B,  C,  D)  and  the  wave  or  the  sediment 
characteristics  of  a beach  profile.  The  best  fit  was  obtained  using  the  dimensionless  fall 
velocity,  £3  = //,/ wT - For  the  range  of  l .5  < £3  < 5,  the  following  expressions  were 
obtained 

-4  = 0.21-0.02£J  (1.26) 

B = 0.89exp(-l  ,24Q)  (1.27) 

C = 0.66  + 0.04£2  (1.28) 

D = 0.22exp(-0.83fi)  (1.29) 

1. 2.2.3  Non-uniform  grain  size  effects 

On  most  beaches,  the  median  grain  size  decreases  with  distance  offshore, 
producing  a profile  that  is  steeper  inshore  and  more  gently  sloping  offshore.  Based  on 
the  same  assumption  as  Dean's  equilibrium  profile  concept,  a modified  equilibrium 
profile  model  was  proposed  for  varying  sediment  size  along  the  beach  by  Larson  (1991). 
A simple  description  of  the  variation  in  equilibrium  energy  dissipation,  D^,  with  distance 
from  the  shoreline  is  assumed  as 

O,=D.+(D0-D.)cxp(-2y)  (1.30) 

where  Do  and  D«  are  the  equilibrium  energy  dissipation  for  the  material  at  the  shoreline 
and  far  offshore  respectively  (Do  > D,),  and  X is  an  empirical  coefficient  describing  the 
rale  at  which  D„  approaches  D,  . The  equilibrium  profile  shape  is  given  by 


(1.31) 


A=^+i<fL'lxl_e'*))! 

in  which  A is  defined  based  on  D«.  Three  different  beach  profile  data  sets  were  used  to 
compare  the  modified  equilibrium  profile  shape  with  the  classical  2/3-powcr  equilibrium 
profile,  and  a better  fit  with  the  measured  profile  was  found  by  the  modified  equilibrium 
profile  form. 

1.2.3  Forces  Acting  in  the  Nearshore 

An  equilibrium  beach  profile  represents  a balance  of  "destructive"  and 
"constructive"  forces  acting  on  the  beach.  At  present,  neither  the  complete  identification 
nor  the  quantification  of  these  forces  is  well  understood.  Dean  (199S)and  Dean  and 
Dalrymple  (2002)  described  these  forces  in  detail. 

Destructive  forces  arc  those  that  cause  offshore  transport  and  beach  erosion. 
Gravity  is  the  most  obvious  destructive  force.  It  acts  downslopc  reducing  the  profile 
slope.  Gravity  also  tends  to  smooth  any  irregularities  that  occur  in  the  profile.  Another 
significant  destructive  force  is  the  turbulence  generated  by  wave  breaking  within  the  surf 
zone.  The  turbulent  fluctuations  dislodge  sediment  particles  and  suspend  them  into  the 
water  column.  The  seaward  return  flow,  undertow,  then  transports  these  suspended 
sediment  offshore.  The  undertow  also  causes  seaward-directed  shear  stresses  that 
contribute  to  an  offshore  bedload  transport. 

Constructive  forces,  on  the  other  hand,  cause  onshore  transport  and  beach 
accretion.  Three  individual  constructive  forces  are  identified.  The  first  is  due  to  the  net 
onshore  shear  stresses  at  the  bottom  resulting  from  the  nonlinear  shallow-water  wave 
forms.  The  wave  crests  arc  higher  and  of  less  duration  than  the  troughs  for  the  nonlinear 
waves.  The  second  constructive  force  is  that  due  to  the  streaming  velocities  in  the 


direction  of  propagating  waves  within  the  bottom  boundary  layer.  An  additional 


constructive  force  is  due  to  the  intermittent  suspension  and  selective  transport  of  the 
particles  by  the  wave-induced  crest  velocities.  This  force  could  be  constructive  or 
destructive  depending  on  the  transport  direction  caused. 

1.2.4  Parameters  Affecting  Profile  Evolution 
1.2.4.1  Sand  grain  size 

The  transportability  of  a sand  grain  is  affected  by  its  characteristics  (such  as  size, 
shape,  porosity,  specific  gravity,  etc.).  A hydraulically  relevant  sand  particle  parameter 
its  fall  velocity.  Coarse  sand  grains  (which  fall  rapidly)  are  less  likely  than  finer  sand 
grains  to  be  suspended  in  the  water  column.  Coarse  sand  also  can  sustain  larger  shear 
stresses  without  moving  than  fine  sand.  Thus,  coarse  sediments  can  withstand  the 
destructive  forces  more  easily  than  fine  sands,  and  usually  form  a steeper  beach  slope. 
Figure  1-5  compares  different  sediment  size  effects  on  profile  shape  for  the  same  wave 


1. 2.4.2  Wave  height 


For  the  same  profile,  increasing  wave  height  results  in  a wider  surf  zone  and  more 
wave  energy  dissipation,  thus  causing  more  erosion.  Under  higher  waves  such  as  storm 
conditions,  a so-called  “storm"  or  "winter"  profile  will  be  developed;  single  or  multiple 
bars  usually  will  appear  on  this  kind  of  profile.  By  contrast,  under  lower  waves,  a 
"normal"  or  "summer"  profile  will  be  formed.  Transport  in  these  two  conditions  appears 
to  occur  in  significantly  distinct  modes  and  with  markedly  disparate  time  scales;  as  a 
result,  the  difficulties  in  predictive  capabilities  differ  substantially.  Figure  1-6  compares 
different  wave  height  effects  on  the  profile  evolution. 


Figure  1-6.  Comparison  of  profile  evolution  under  different  wave  heights  (Data  from 
Kraus  and  Larson,  19S8). 

1.2. 4.3  Wave  period 

According  to  the  fall  time  model  of  Dean  (1973),  if  the  fall  lime  of  a suspended 
sand  particle  is  less  than  half  the  wave  period,  the  sand  grain  should  move  onshore. 
Alternatively,  if  the  fall  time  is  greater  than  half  the  wave  period,  the  sand  particle  would 
be  carried  offshore.  Thus  longer  wave  periods  tend  to  cause  sediment  to  be  transported 


shoreward  and  the  mean  shoreline  to  advance  seaward  such  that  the  average  beach  slope 


is  steeper  than  for  shorter  period  waves.  Figure  1 -7  shows  different  profile  development 
under  different  wave  periods.  In  both  tests,  grain  sizes  are  the  same,  and  the  incident 
wave  heights  are  almost  the  same.  The  test  with  longer  wave  period  has  an  accretionary 
profile;  on  the  other  hand,  the  lest  with  shorter  wave  period  has  an  erosional  profile. 


1 .2.4.4  Regular  and  irregular  waves 

Random  waves  typically  break  across  a much  wider  portion  of  the  beach  profile 
than  monochromatic  waves  with  the  same  average  energy,  producing  a more  evenly 
distributed  energy  dissipation  across  shore.  Thus,  the  net  cross-shore  transport  rate 
distribution  is  expected  to  be  more  uniform  under  random  waves,  causing  a more  gradual 
evolution.  Figures  1 -8  and  1 -9  show  profile  changes  and  transport  rates  of  two  runs  in 
SUPERTANK  (Kraus  and  Smith,  1994)  under  monochromatic  and  random  waves.  The 
beach  profile  was  first  exposed  to  random  waves  until  near  equilibrium.  Then  a 


monochromatic  wave  was  used,  with  a wave  height  the  same  as  the  significant 1 


height  of  the  random  wave  and  a wave  period  the  same  as  the  peak  spectral  period.  The 
shape  of  the  bar  clearly  changed  after  it  was  exposed  to  monochromatic  waves.  Bar  size 
increased  simultaneously  as  a pronounced  trough  developed. 


Figure  1-9.  Transport  rates  under  random  and  monochromatic  waves  for  case  ST-10  in 
SUPERTANK  (H™  = 0.8  m,  Tp  = 3 s). 


1.3  Objectives  and  Scope  of  Study 

At  present,  although  progress  has  been  made  in  cross-shore  sediment  transport 
modeling,  there  are  no  known  models  which  have  been  demonstrated  to  predict 
adequately  the  time  scales  for  nourishment  profile  equilibration  for  practical  purposes. 
The  main  objective  of  this  study  is  to  provide  a better  understanding  of  cross-shore 
sediment  transport  relationships  and  develop  a numerical  model  to  represent  the  beach 
profile  equilibration  process  and  determine  the  time  scale  of  profile  evolution,  which  can 
be  applied  by  practicing  engineers. 

Chapter  2 provides  a detailed  literature  review  on  the  subject,  which  includes 
previous  valuable  findings  in  laboratory  experiments,  and  some  analytical  and  numerical 
models  developed  previously. 

Chapter  3 is  the  main  pan  of  the  dissertation  and  presents  detailed  descriptions  of 
the  model  developed  herein.  Three  different  wave  energy  dissipation  models  under 
monochromatic  waves  arc  examined:  (I)  H = jh , the  spilling  breaker  assumption  which 
indicates  a saturated  surf  zone,  (2)  the  wave  transformation  model  of  Dally,  Dean  and 
Dalrymple  (1985),  and  (3)  Dally,  et  al  model  coupled  with  the  roller  model  of  Dally  and 
Brown  (1995).  The  relationship  between  the  transport  rales  and  the  deviation  of  the 
actual  wave  energy  dissipation  per  unit  volume  from  equilibration  is  adopted.  An 
empirical  relationship  is  developed  which  incorporates  wave  effects  in  the  profile  scale 
parameter,  A (see  Eq.1.1),  of  the  equilibrium  beach  profile.  An  implicit  finite  difference 
scheme  with  a uniformly  spaced  grid  in  the  offshore  direction  is  applied  to  solve  the 
model  governing  equations.  The  transport  and  profile  evolution  under  random  waves  are 
determined  using  a wave-by-wave  approach. 


The  model  calibrations,  evaluations  and  sensitivity  analyses  arc  given  in  Chapters  4 
and  5.  The  data  from  lour  large  wave  flume  experiments  are  utilized  for  the  model 
calibration  and  evaluation,  these  data  include  both  regular  waves  and  random  waves. 
Sensitivity  analyses  investigate  the  influences  of  several  model  parameters  on  the  model 
performance,  [n  Chapter  6,  both  additional  laboratory  data  and  field  data  arc  applied  to 
further  investigate  the  model  performance.  The  model  comparisons  with  SBHACH 
(Larson  and  Kraus,  1 989)  and  CROSS  (Zheng,  1996;  Zheng  and  Dean,  1997)  are  also 
provided  in  Chapter  6. 

Chapter  7 presents  the  summary  and  conclusions  from  the  study  and  proposes  some 
recommendations  for  future  model  development.  The  program  code  in  FORTRAN90  is 
presented  in  the  appendix,  where  more  case  studies  are  also  provided. 


CHAPTER  2 
LITERATURE  REVIEW 

2.1  Laboratory  Studies 

Due  to  the  complexity  of  forcing  and  beach  change  processes  in  nature,  laboratory 
examination  under  controlled  conditions  remains  a powerful  tool.  Many  laboratory 
research  studies  have  been  conducted  during  past  decades,  including  small-scale  and 
large-scale  wave  flume  experiments.  These  experiments  provide  useful  information  and 
guidance  for  our  understanding  and  modeling  of  beach  proflle  evolution  processes. 

Rector  (1954)  conducted  laboratory  experiments  of  equilibrium  beach  profile 
characteristics.  He  suggested  that  the  wave  steepness  and  the  ratio  of  median  grain  size 
to  wave  length  (djo/Lo)  played  an  important  role  in  beach  profile  formation. 

Walts  (1954)  investigated  the  influence  of  varying  wave  periods  on  beach  profiles. 
It  was  found  that  the  final  foreshore  and  offshore  slopes  produced  were  essentially 
identical  under  constant  period  wave  and  varying  period  waves  (with  the  mean  value 
equal  to  the  constant  period  wave);  and  tile  variable  period  tests  created  less  actual 
material  movement  and  less  apparent  bed  features  (bar  and  trough)  than  the  constant 
period  tests. 

Raman  and  Earattupuzha  (1972)  examined  the  stable  points  on  beach  profiles 
commencing  with  a planar  profile.  A stable  point  is  the  point  on  a beach  profile  that  does 
not  undergo  any  considerable  change  during  the  process  of  profile  modification  due  to 
incident  waves  of  essentially  constant  characteristics.  Figure  2-1  shows  examples  of 
stable  points  on  storm  and  normal  profiles.  For  a storm  profile,  there  exists  a stable  point 
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near  die  breaker  zone  and  a second  stable  point  may  exist  in  the  offshore  zone  depending 
on  the  steepness  of  the  incident  waves.  In  the  case  of  a normal  profile,  the  stable  point 
does  not  exist  in  the  strict  sense,  but  there  is  one  point  on  the  modified  profile  which 
seems  to  serve  the  same  function  as  a stable  point  in  the  storm  profile  case.  They  found 
that  the  depth  at  the  stable  point  was  proportional  to  the  dimensionless  parameter  HogT/v, 
where  v is  the  kincmatical  viscosity  of  water. 


•^=C(tan/»)-‘’J,(%)0 


(2.1) 


where  p is  the  initial  slope.  For  Type  I profile:  C a 4,  Type  n profile:  4£C  <S,  Type 
m profile:  C£8. 

Watanabe  et  al.  (1980)  conducted  a series  of  small-scale  experiments  using  two 
different  sand  diameters  (0.2  mm  and  0.7  mm).  Each  experiment  began  with  a uniform 
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initial  slope  of  1/10  or  1/20  and  was  conducted  under  different  wave  conditions  for  one- 

hour.  They  confirmed  the  beach  profile  classification  criterion  of  Sunamura  and 

Horikawa  (1974)  and  related  the  net  sediment  transport  rates  to  the  Shields  parameter  by 

® = 4(4<-4'„)  (2.2) 

where  d>  = is  the  dimensionless  transport  rate,  w the  fall  velocity,  d the  grain 

diameter,  4*  = is  the  Shields  parameter  and  '!'  is  the  critical  value  for  the 

p(s-\)g<i 

initial  sediment  movement.  Figure  2-3  shows  the  results  of  their  relationship  between  net 
transport  rate  and  Shields  parameter. 


Figure  2-2.  Three  types  of  beach  profiles  (Sunamura  and  Horikawa,  1974) 
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Kajima  el  al.  (I9S2)  classified  the  nel  transport  rate  into  five  types  (Figure  2-4) 
based  on  their  large  wave  flume  data.  They  also  proposed  a model  of  beach  profile 
change  by  assuming  that  the  total  on-oflshorc  transport  rate  Q could  be  represented  as  the 
sum  of  three  independent  components  of  net  sand  transport  rate 

e-a-a+a  (2.3) 

where  Qi  = net  onshore  transport  outside  the  surf  zone  caused  by  nonlinearity  of  the  fluid 
motion  under  waves;  Q2  = net  offshore  transport  inside  the  surf  zone  caused  by  time- 
averaged  seaward-directed  currents;  and  Q;  = nel  onshore  transport  in  the  swash  zone 
caused  by  wave  run-up.  The  three  components  were  each  assumed  to  be  Gaussian 
distributed. 


Figure  2-4.  Types  of  net  transport  rate  distributions  (Kajima  et  al.,  1982). 

Sunamura  and  Maruyama  (1987)  used  large  wave  flume  experimental  data  to 
examine  eroding  beach  characteristics.  They  found  the  bar  migration  directions  could  be 
classified  by  dimensionless  parameter  K defined  as 


lfK>  10,  bars  move  offshore;  whereas  they  move  onshore  if  K<  10. 


The  lime  required  for  the  attainment  of  steady  slate,  t,,  in  the  height  of  seawatd 


migrating  bars,  I,,  was  found  to  be  described  by 


(2.5) 


where  the  parameter  in  the  bracket  is  the  square  of  the  Iribarren  number  and  indicates  the 
breaker  type.  With  increasing  values  of  this  parameter,  the  breaker  type  changes  from 
plunging  to  spilling.  The  time  variation  of  the  bar  crest  location,  X,  was  expressed  by  an 
exponential  form  as 


where  Xc  is  the  equilibrium  value  of  X,  T is  a so-called  equilibrium  time  scale.  At  t=r , 
X reaches  63%  of  the  total  change  Xc.  T was  found  to  be 


Other  experiments  include  the  large  wave  tank  (LWT)  experiment  in  1956-1957 
and  1962  by  Saville  in  Dalecarlia,  Washington,  D.C.  (Kraus  and  Larson,  1988),  Large 
Wave  Flume  in  Hannover,  Germany  (Dettc  and  Uliczka,  1987;  Dettc  ct  al„  1998), 
experiments  at  the  Central  Research  Institute  of  Electric  Power  Industry  in  Japan 
(Sunamura  and  Maruyama.  1987),  SUPERTANK  project  in  the  large  wave  tank  at  the 
Otto  Hinsdale  Wave  Research  Laboratory,  Oregon  State  University  (Kraus  and  Smith, 

1 994),  and  LIP  1 1 D Delta  Flume  experiments  (Roevink  and  Renicrs,  1 995),  etc.  Some  of 
the  tests  include  regular  waves,  whereas  others  include  both  regular  and  irregular  waves. 
They  provide  many  useful  data  for  later  studies  on  the  profile  evolution  process. 


* = Ar,(l-exp  (~/T)) 


(2.6) 


r = 2.87xl0’ 


(2.7) 
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2.2  Numerical  Modeling  Studies 

Some  recent  interesting  reviews  of  cross-shore  sediment  transport  models  have 
been  reported  from  different  methods  and  viewpoints  by  Broker  Hedegaard  et  ol.  (1992), 
Stcetzel  (1993),  and  Schoonees  and  Theron  (1995).  Broker  Hedegaard  et  al.  made  an 
inlercomparison  of  six  different  coastal  profile  models  against  measured  profile 
evolutions  from  a large  wave  flume  (under  dune  erosion  conditions).  Stcetzel  classified 
transport  formulae  as  empirical  transport,  shear  stress  related  transport,  energy  dissipation 
related  transport,  energy  dissipation  related  suspension  and  velocity  related  transport. 
Schoonees  and  Theron  summarized  ten  internationally  developed  transport  models 
according  to  their  theoretical  basis  and  morphological  verification  data. 

Generally,  numerical  models  of  profile  evolution  can  be  termed  as  "open  loop"  and 
“closed  loop"  (Dean,  1995;  Dean  and  Dalrymple,  2002).  An  open  loop  model  is  not 
constrained  to  a final  profile  and  the  sediment  transport  is  determined  by  sediment  fluid 
properties.  Open  loop  models  required  detailed  or  idealized  representation  of  the 
hydrodynamics  and  sediment  transport  processes.  Closed  loop  models  are  characterized 
by  the  specification  of  a particular  "target  profile,"  to  which  the  computed  profile  will 
converge  if  the  forcing  mechanisms  are  held  constant  indefinitely.  Cross-shore  transport 
is  modeled  by  deviations  of  a beach  profile  from  equilibrium.  Open  loop  models  arc 
microscale  related  approaches,  whereas  closed  loop  models  are  more  macroscale  and  use 
empirical  equations  based  on  laboratory  and/or  field  data  to  provide  reasonable  solutions. 

2.2.1  Closed  Loop  Models 

2.2.1. 1 EDUNE  model 

One  of  the  first  closed  loop  models  developed  was  that  of  Kriebel  (1982)  and 
Kricbcl  and  Dean  (1985)  (originally  called  EBEACH  and  later  called  EDUNE)  based  on 


the  energy  dissipation  concept.  The  basic  concept  was  that,  for  a uniform  sand  size 


across  the  profile,  a uniform  energy  dissipation  rate  per  unit  volume  resulted  in  the  Ay' 
profile  as  an  equilibrium  beach.  If  the  beach  profile  differed  from  this  equilibrium  form, 
the  dissipation  rates  across  the  surf  zone  would  differ  from  the  equilibrium  value.  The 
cross-shore  sediment  transport  rate  at  any  point  in  the  surf  zone,  Q,  was  related  to  the 
difference  between  the  actual  and  equilibrium  levels  of  wave  energy  dissipation  in  the 
surf  zone  (D  and  D*,  respectively). 

Q = K{D-D.)  (2.8) 

where  K is  a transport  rate  parameter  found  during  calibration  tests  first  based  on  the 
LWT  data  of  Saville  and  later  Hurricane  Eloisc  profile  data  (Kriebel,  1986),  resulting  in  a 
value  ofS.TSxlO-4  m4/N.  Using  shallow-water  wave  theory  and  assuming  spilling 
breakers,  D was  determined  as 


(2.9) 


where  y = 0.78  is  the  breaker  index,  p is  density  of  sea  water,  and  g is  gravity.  To  obtain 
time-dependent  profile  evolution,  an  implicit  double-sweep  method  was  employed  to 
solve  the  transport  equation  and  the  conservation  of  sand  equation: 


dy_  SQ 
ill  dh 


(2.10) 


Both  Q and  y are  defined  as  positive  in  the  offshore  direction  and  h is  positive  below  the 
still  water  level.  OITshore,  sediment  transport  was  assumed  to  be  effectively  limited  at 
the  break-point  of  the  incoming  waves.  Later  developments  included  the  effects  of  wave 
run-up,  dune  steepening,  and  dune  overwash  in  the  onshore  boundary  region  (Kriebel, 


1990).  Fig 


i 2-5  and  2-6 ! 


aplcs  of  applic 


for  large 


lank  and  field  data,  respectively. 


2.2.1. 2 SBEACH  model 


Based  on  observations  from  large  wave  tank  data  (mainly  from  Kraus  and  Larson, 
1988  and  Kajima  et  al„  1982),  Larson  (1988)  and  Laison  and  Kraus  (1989)  developed  a 
model  called  SBEACH  to  predict  beach  profile  change  resulting  from  cross-shore  sand 
transport,  focusing  on  the  main  morphologic  features  of  bars  and  berms.  Changes  in  the 
beach  profile  were  assumed  to  be  produced  by  breaking  waves;  therefore,  the  cross-shore 


lined  from  I 


properties,  and  the  conservation  equation  of  sand  in  the  cross-shore  direction  was  solved 
to  compute  profile  changes  as  a function  of  time. 


This  model  included  a wave  module  to  calculate  the  wave  height  and  other  wave- 
related  properties  at  grid  points  on  a profile  line.  The  wave  height  was  determined  by 
linear  wave  theory  in  regions  of  nonbreaking  by  assuming  locally  plane  and  parallel 
depth  contours,  and  using  the  Dally,  Dean  and  Dalrymple  ( 1 984)  model  in  the  surf  zone. 
Based  on  the  wave  dissipation  model  ofDallyetal.,  an  equilibrium  profile  (Eq.2.1 1)  was 


where  m is  shoreline  slope,  the  A.  value  is  different  from  the  A value  used  in  the  usual 
equilibrium  profile  equation  due  to  the  additional  linear  slope  dependence. 

The  distribution  of  transport  rate  across  the  nearshore  zone  was  calculated  using 
different  relationships  in  four  distinct  zones  (Figure  2-7): 

• Zone  I:  Pre-breaking  zone,  extending  offshore  from  the  breaking  point. 

The  sediment  transport  rate  was  expressed  as 


where  Qnp  is  the  sediment  transport  at  the  break  point,  yap  is  the  break  point  location, 
is  die  spatial  decay  coefficient  given  by 


(2-11) 


0=e«pe>tp[-2iO'-p«p)] 


(2.12) 


where  d)o  is  the  median  grain  size  in  mm,  Hb  is  the  breaking  wave  height  in  m,  and  the 
unit  of  ki  ism*1. 


—if 
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where  Ho  and  Lo  are  ihc  deep  walcr  wave  height  and  length  respectively,  w is  the 
sediment  fall  velocity  and  T is  the  wave  period. 


Figure  2-7.  Principal  zones  of  cross-shore  transport  (Larson  and  Kraus,  1 989). 

Later,  Wise  ct  al.  (1996)  modified  the  model  to  simulate  cross-shore  profile 
changes  under  random  waves.  A semi-analytic  random  wave  model  based  on  Dally  et  al. 
( 1 984)  model  was  developed  to  predict  the  wave  height,  the  fraction  of  broken,  unbroken 
and  reformed  waves,  and  the  wave  energy  dissipation  across  shore.  Then  the  original 
transport  equations  were  modified  to  represent  transport  under  broken  and  unbroken 
waves.  An  algorithm  for  dune  erosion  produced  by  overwash  was  also  developed  in  this 
version.  Later,  Larson  and  Kraus  (1998)  included  hard  bottom  effects  into  the  model. 
2.2.1 3 CROSS  model 

Zheng  ( 1 996)  and  Zheng  and  Dean  (1997)  developed  a non-linear  cross-shore 
model  (CROSS),  which  was  a modification  of  the  EDUNE  model  and  based  on  scale 
analysis. 

From  Eq.2.8,  the  scaling  relationship  in  Eq.2.17  can  be  established. 


(2.17) 
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in  which  the  subscript  r indicates  a ratio  of  model  to  prototype  values.  For  an  undistorted 
model,  according  to  the  definition  of  D (Eq.2.9),  the  disequilibrium  scale  ratio  can  be 
expressed  in  terms  of  the  length  scale  ratio,  Lr,  as 

(D-D.\=jr,  (2.18) 

On  the  other  hand,  the  Froudc  relationship  yields  the  lime  scale  ratio,  Tr , as 

(2.19) 

From  the  relationships  presented,  the  cross-shore  sediment  transport,  Q,  should  be  scaled 

Q.=y  = L?‘  (2-20) 

If  K is  considered  to  be  a constant,  the  scaling  ratio  of  the  Kriebel  and  Dean  formulation 
can  be  shown  to  be 

a =V*  (2.2i) 

which  does  not  provide  a valid  scaling  of  the  transport  unless  K varies  with  scale. 
A transport  rate  equation  was  then  proposed 

Q*‘K(D-D.p-D.\"'  (2.22) 

which  results  in  the  following  scale  relationship 

Q,  = K,(D-  D.  ),|D  - D.|;-  = K,L-n  (2.23) 

n = 3 was  determined  by  assuming  K,  independent  of  the  length  scale  and  equal  unity 
such  that  both  the  scaling  relationship  and  convergence  to  the  equilibrium  profile  are 
satisfied.  They  argued  that  this  relationship  could  provide  a reasonable  explanation  for 


I 

i 
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Figure  2-9.  Comparison  of  linear  (n=l)  and  nonlinear  (n=3)  transport  rates  with  measured 
data  in  German  large  wave  flume  (Zheng,  1996). 

2.2.1.4  Kobayashi  model 

Kobayashi  (1987)  suggested  a two-parameter  sediment  transport  rate  equation  as 
following 

g = D(A)^-A';  OsAsA,  (2.24) 

where  D and  K arc  empirical  functions  and  related  to  the  wave  condition  and  sediment 
characteristics  inside  the  surf  zone.  The  suggested  form  of  D and  K by  Kobayashi  (1987) 

D(h) =ah''~‘v‘  ,K  = naA'"  (2.25) 

which  generate  a steady  stable  solution,  A = A(x  + R)" , similar  to  the  equilibrium  beach 
profile  by  Dean  (1977).  A and  a are  parameters  related  to  the  sediment  size  and  fall 
velocity.  When  n =2/3,  Eq.2.24  becomes  identical  to  that  ofKriebel  and  Dean  (1985). 

Lee  et  al.  (1996)  used  this  formulation  and  developed  a numerical  model  to 
simulate  the  cross-shore  profile  evolution  with  moving  boundaries.  Further,  Lee  ct  al. 
(1999)  extended  this  model  to  a multi-domain  hybrid  numerical  model  including  the 
generation  and  movement  of  nearshore  bars  by  breaking  waves.  Similar  to  the  SBEACH 
model  by  Larson  and  Kraus  (1989),  the  computation  domain  is  divided  into  three 
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subdomains,  including  post-breaking  zone,  breaking  zone  and  pre-breaking  zone.  In  each 
region,  different  sediment  transport  rate  equations  were  used. 

In  the  post-breaking  zone, 

Q - D~  K , - m £x  £ 1,(0  (2.26) 

D = aJh^,K  = laA''!  (2,27) 

where  hn  and  xp  are  the  water  depth  and  distance  from  the  origin  to  the  plunge  point,  R(t) 
is  the  distance  from  the  origin  to  the  point  h(x,t)  = 0.  In  the  breaking  zone  and  pre- 
breaking  zone,  the  sediment  transport  rates  were  assumed  to  have  the  following  explicit 
expressions: 

Q,  = r/.(r)cxp(-2,(x-x.)'),  x„  £x£xt  (2.28) 

Q,  = ?, (f ) exp(-2, (x - ,ts )) , xt  £x£*>  (2.29) 

where  q.(l)  and  </,(f)  are  time-dependent  amplitude  functions  for  the  respective 
sediment  transport  rates  in  breaking  zone  and  pre-breaking  zone.  Eq.  2,28  is  different 
from  the  exponential  function  used  in  SBEACH.  This  equation  means  that  the  sediment 
transport  rate  in  breaking  zone  has  a Gaussian  distribution  (Kajima  et  al.  1982)  which  has 
a maximum  value  at  a point  x-  positioned  between  plunging  and  breaking  points.  The 
relationship  between  /,  and  2,  was  found  to  be 

*,=£  (2.30) 

The  empirical  formula  for  2,  in  the  model  was 


(2.31) 
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The  model  was  formed  by  the  combinalion  of  numerical  solutions  (post-breaking 
zone  I)  and  empirically  based  analytic  solutions  (breaking  zone  II  and  pre-breaking  zone 
III).  This  model  can  straightforwardly  be  extended  to  multi-bar  problems  by 
systematically  adding  more  subdomains.  Figures  2-10  and  2-1 1 show  comparisons  of 
single-bar  and  multi-bar  results  with  LWT  experimental  data  (Kraus  and  Larson,  1988). 


Q(m'/m/hr) 


Distance  offshore(m) 


Figure  2-10.  Comparison  of  the  averaged  net  sediment  transport  rates  computed  from  the 
single-bar  and  multi-bar  numerical  models  with  CE  Case  300  experimental 
data  after  5 hours  (measured  profile  had  two  bars)  (Lee  et  al.,  1999). 


Two  analytical  or  semi-analytical  closed  loop  models  have  been  documented  in  the 
literature.  Kobayashi  (1987)  combined  the  conservation  equation  and  the  transport 
relationship  in  Eq.2.24  into  one  equation  representing  the  profile  response. 


(2.32) 
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where  D and  K are  defined  in  Eq.2.25.  The  landward  and  seaward  limits  of  ihe  domain 
were  treated  as  moving  boundaries  and  the  eroding  dune  face  was  represented  as  a 
vertical  face.  An  analytical  solution  was  determined  for  a particular  initial  condition.  It 
was  found  that  the  problem  is  mathematically  similar  to  a difiusion  problem  in  which  the 
convergence  to  the  steady  state  problem  first  occurs  rapidly  and  later  approaches  the 
steady  state  solution  asymptotically.  The  model  results  were  found  to  approximate  the 
numerical  solutions  of  the  EDUNE  model. 

Kricbel  et  al.  (1991)  and  Kricbcl  and  Dean  (1993)  developed  and  illustrated  a 
convolution  method  of  representing  beach  profile  evolution.  The  method  is  based  on  the 
observation  from  numerical  solutions  that  for  a constant  forcing  function,  i.e.,  water 
level,  a particular  contour  y(t)  tends  to  approach  the  equilibrium  value  y,  asymptotically 
in  approximate  accordance  with  the  following  relationship 

J'(')  = ^.(l-e‘"';)  (2.33) 

in  which  7)  is  the  characteristic  time  scale  of  the  exponential  response.  The  solution  is 
obtained  in  the  form  of  a convolution  integral  of  the  response  function  and  the  water  level 
history,  fit), 

3'W  = YLJ/(r)exp[-(r-r)/ri]<fr  (2.34) 

The  method  was  found  to  yield  very  similar  results  to  those  obtained  through  numerical 
modeling. 


2.2.2  Open  Loop  Models 


Open  loop  models  usually  relate  cross-shore  sediment  transport  to  the  detailed 
physics  of  the  flow  field  such  as  sediment  concentration,  fluid  velocity,  and  bottom  shear 
stress,  etc.  Some  of  these  models  will  be  reviewed  briefly  here. 

2.2.2. 1 Dally  and  Dean  model 

The  Dally  and  Dean  model  (Dally,  1980;  Dally  and  Dean,  1985)  was  based  on  net 
time-averaged  flux  of  suspended  sediment  past  a section  in  the  nearshore  zone  given  by 


where  u(z)  = average  horizontal  velocity  at  level  z,  and  C(z)  = suspended  sediment 
concentration  at  level  z. 

The  flow  regime  was  divided  into  an  upper  layer  where  only  mean  flow  was 
considered  and  a lower  layer  where  both  mean  flow  and  orbital  velocity  were  taken  into 
account.  The  interface  between  layers  was  determined  by  the  elevation  at  which  the 
assumed  uniform  sediment  would  fall  back  to  the  bottom  in  one  wave  period.  Linear 
wave  theory  was  used  to  predict  orbital  velocities  and  stream  function  wave  theory  had 
been  applied  to  obtain  second-order  mean  flow  velocities.  The  cross-shore  wave  heights 
were  obtained  from  die  breaker  decay  model  of  Dally  ct  al.  ( 1 984). 


The  suspended  sediment  concentration  was  considered  to  be  exponential  as  given 
by 


where  C4  is  the  concentration  at  an  arbitrary  reference  elevation,  : , . The  ratio  of  the 


Q = Ju(r)C(z)<& 


(2.36) 


(2.37) 


sediment  fall  velocity,  w, 


i , governed  the  exponential 
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nature  of  the  profile.  The  shear  velocity  was  assumed  to  be  the  sum  of  the  shear 
velocities  due  to  wave-induced  bottom  shear  and  breaking-induced  turbulence. 

The  continuity  equation  was  then  solved  in  an  explicit  finite  difference  scheme. 
However,  to  prevent  numerical  instability  just  outside  the  surf  zone,  empirical  transport 
spreading  was  applied.  The  model  could  produce  both  normal  and  storm  profiles  (Figure 
2-12),  however,  the  beach  profile  shapes  were  not  very  favorable  quantitatively. 

Stcetzcl  ( 1 993)  also  presented  his  cross-shore  transport  model  as  the  integral  of  the 
product  of  concentration  and  velocity  over  time  and  depth.  The  sediment  concentration 
was  found  from  a one-dimensional,  non-stationaiy  convection-diffusion  equation.  The 
reference  concentration  at  the  bed  was  assumed  to  be  related  to  the  breaking  type  and  the 
dissipation  of  the  turbulent  kinetic  energy.  The  cross-shore  transport  rate  was  calculated 
in  two  parts,  that  is,  below  and  above  the  mean  trough  level.  Above  the  trough  level, 
concentration  was  a constant  and  given  by  that  at  the  mean  water  level. 


The  model  included  a special  treatment  of  the  wave  transition  zone.  The  transition 


zone  length.  /, , was  calculated  by  a "length  method". 


/,  = 0.556m,  £,#-'■ 


(2.39) 


Wo" 


and  m„  Lt,  //,  are  respectively  the  bed  slope,  wave  length  and 


wave  height  at  breaking. 

After  the  transition  zone  length  was  given,  an  expression  for  the  fraction  of  waves. 
Q, . in  the  breaking-induced  turbulence  region  was  determined  as 


where  Qh  is  the  fraction  of  wave  breaking  in  the  surf  zone  predicted  by  a wave 
transformation  model,  y is  the  cross-shore  coordinate.  In  the  Naim  and  Southgate  model, 
fixed  surface  and  credible  surface  problems  were  also  addressed.  The  model  has  been 
compared  with  both  small  scale  and  prototype  laboratory  data,  and  the  rate  and  type  (i.e. 
erosion,  deposition  and  bar  migration)  of  beach  response  can  be  well  estimated. 

Plant  el  al.  (200 1 ) developed  a simple  parametric  model  for  the  cross-shore 
transport  based  on  Bagnold's  model. 


where  c,  r, , and  r3  arc  free  parameters  to  be  calibrated,  H is  the  local  wave  height.  The 
first  torm  in  parentheses  is  a wave  stirring  term,  and  the  terms  in  braces  form  a 
rectification  (unction  representing  the  influences  of  downslope  transport  onshore 
transport  due  to  weak  wave  nonlinearity  (e.g,  wave  steepness),  and  offshore  transport  due 
to  strong  nonlinearity  (e-g.  undertow). 


aw=a(p-/,) 


(2.40) 


(2.41) 


2.2.2.3  Watanabc  model 


Walanabc  (1988)  preaenled  a sediment  transport  rate  equation  related  to  local 
wave-current  conditions  and  bottom  shear  stress.  The  total  transport  rate  was  the  sum  of 
contributions  from  mean  current,  Qr , and  waves,  Q„,  which  are  expressed,  respectively, 
by 

Q,  = A,(r -r,p  I pg  (2.42) 

Q.~A.F„(t-x,)umlpg  (2.43) 

where  A,  and  A,  are  non-dimensional  coefficients,  U the  current  velocity  vector,  um  the 
maximum  near-bottom  orbital  velocity  vector,  F„  a direction  function  for  wave-induced 
net  transport,  r the  maximum  bottom  shear  stress  in  a wave-current  coexisting  system, 
and  r,  the  critical  shear  stress  for  the  onset  of  the  general  movement. 

These  formulas  are  based  on  the  wave  power  concept  and  assume  that  the 
sediments  set  in  motion  by  the  excess  shear  stress  under  combined  wave-current  action, 
arc  transported  by  both  mean  currents  and  wave  motion  into  the  respective  directions. 
The  coastal  profile  model  WATAN3  developed  by  University  of  Liverpool  (reviewed  in 
Broker  Hcdegaard  et  al„  1992)  is  based  on  the  Watanabe  transport  equations. 

There  are  some  other  transport  relationships  such  as  the  one  developed  by 
Nishimura  and  Sunamura  (1986),  in  which  the  transport  rate  was  related  empirically  to 
the  Urscll  number  and  the  Hallermeicr  parameter,  which  represent  the  skewness  of  water 
particle  velocity  profile  and  the  intensity  of  sediment  movement,  respectively.  Camcncn 
and  Larroude  (2003)  compared  four  process-related  sediment  transport  formulas  with 
field  and  experimental  data  and  gave  some  recommendations  on  the  use  of  these  different 


formulas. 


Most  open  loop  profile  change  models  include  modeling  of  hydrodynamic 
processes  to  different  degrees.  The  sediment  transport  calculations  are  related  to  the 
properties  of  sediment  and  flow  field.  Van  Rijn  ct  al.  (2003)  compared  five  major  coastal 
profile  change  models  developed  in  Europe.  The  five  models  are:  UNIBEST-TC  of  Delft 
Hydraulics.  COSMOS  of  HR  Wallingford,  CROSMOR2000ofUnivetsityofUtrecht, 
BEACH1/3D  of  University  of  Liverpool,  and  CtlRC  of  University  of  Catalunya.  The 
computed  results  from  these  models  have  been  compared  with  large  wave  flume  data  and 
field  data.  Figures  2-13  and  2-14  show  the  calculated  profile  evolution  and 


hydrodynamic  results,  respectively,  with  comparison  to  the  measurements  in  one  case 


from  COAST3D-Storm  Experiment  at  3D  field  site  of  Egmond,  Netherlands.  It  is  noted 
that  the  all  five  models  present  good  fits  for  the  profile  although  the  calculated 
hydrodynamics  are  quite  different.  It  seems  that  considerable  uncertainty  in  the 
understanding  of  surf  zone  hydrodynamics  and  sediment  transport  processes  still  exists. 


Figure  2-14.  Computed  significant  wave  heights,  longshore  current  velocities,  cross-shore 
current  velocities  and  peak  orbital  velocities  by  UNIBEST-TC,  COSMOS, 
CROSMOR2000,  BEACHI/3D,  and  CURC  with  comparison  to  measured 
results  (Van  Rijn  ct  a!.,  2003). 


and  storm  profiles  can  be  generated  by  most  of  the  open  loop  models.  However,  open 
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stable  for  long  run  limes.  Due  to  the  considerable  complexities  and  uncertainties  of 
hydrodynamics  and  related  sediment  transport  processes  in  the  nearshore  region,  and  the 
uncertainties  in  convergence,  it  appears  that  the  open  loop  models  should  be  considered 
in  the  research  rather  than  applications  arena. 

The  transport  equation  for  the  closed  loop  models  is  chosen  to  ensure  equilibration 
for  fixed  forcing  conditions.  These  closed-loop  models  provide  limited  successful  and 
encouraging  results  and  easy  application  to  practical  engineering  problems.  These 
models  arc  relatively  simple  and  do  not  require  detailed  calculation  of  hydrodynamics 
and  sediment  transport  processes  as  in  open-loop  models.  They  are  all  based  on 
equilibrium  beach  profile  concepts  and  capture  one  of  the  dominant  causes  of  sediment 
transport  inside  the  surf  zone,  namely  wave  energy  dissipation. 

However,  there  are  some  limitations  to  these  closed  loop  models.  (I)  The  grid 
system  used  is  different.  EDUNE  and  CROSS  use  a grid  with  h and  t as  the  independent 
variables,  while  the  SBEACH  and  Lee  ct  al.  models  use  a grid  with  y and  t as  the 
independent  variables.  The  disadvantage  of  the  first  type  of  grid  system  is  that  it  has 
difficulties  in  representing  bar  profiles.  (2)  The  transport  direction  in  EDUNE  and 
CROSS  is  determined  on  the  relative  magnitude  of  wave  energy  dissipation  per  unit 
volume  to  its  equilibrium  value;  SBEACH  uses  an  empirical  relationship  to  determine  the 
transport  direction  along  the  whole  profile,  which  can  not  produce  an  intermediate  profile 
between  fully  erosion  and  accretion  (type  II  profile  in  both  Figures  2-2  and  2-4).  The 
other  severe  shortcoming  of  SBEACH,  as  Dean  (1995)  pointed  out,  is  that  the  transport 
direction  empirical  relationship  used  docs  not  include  any  of  the  profile  characteristics 
and  thus  even  if  the  profile  is  very  steep,  such  as  after  beach  nourishment  and  the  wave 
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and  sediment  characteristics  indicate  ofTshorc  transport,  the  transport  direction  predicted 
by  the  model  is  onshore.  (3)  The  equilibrium  profile  used  in  these  models  depends  only 
on  sediment  size  through  the  profile  scale  factor,  which  cannot  explain  the  seasonal 
profile  slope  changes  in  nature.  (4)  The  wave  energy  dissipation  in  EDUNE  and  CROSS 
is  based  on  the  spilling  breaker  assumption  and  thus  cannot  produce  bar  profiles. 
SBEACH  adopts  a more  robust  wave  model  to  goncratc  bars,  but  deals  with  the  transition 
zone  on  an  empirical  basis.  (5)  The  transport  in  the  swash  zone,  during  overwash  and  in 
front  of  seawalls  are  either  not  included  or  not  well  described  in  previous  models.  (6) 
Another  limitation  is  the  transport  and  profile  changes  under  random  waves.  To  the 
author's  knowledge,  a representative  wave  height  (usually  root-square-mean  wave  height 
or  significant  wave  height)  is  used  in  these  models  except  CROSS  to  simulate  random 

The  motivation  of  this  study  is  to  address  the  limitations  present  in  existing  models 
and  develop  a more  physically  robust  cross-shore  profile  change  model  to  better  represent 
the  cross-shore  sediment  transport  processes  inside  the  surf  zone  (such  as  bar  formation, 
profile  changes  under  random  waves,  etc.),  which  can  be  used  in  practical  field 
applications. 


CHAPTER  3 

CROSS-SHORE  PROFILE  EVOLUTION  MODEL  DEVELOPMENT 
This  chapter  presents  the  cross-shore  profile  evolution  model  developed  in 
conjunction  with  this  dissertation.  The  proposed  model  extends  the  previous  closed  loop 
models  reviewed  in  Chapter  2 and  incorporates  new  findings  to  improve  the  model 
capability  to  better  represent  natural  conditions.  The  following  states  the  questions  to  be 
addressed  in  this  effort. 

• What  is  the  best  representation  of  sediment  transport  rate  in 

Q = K(D-D.)\D-D.\'"  (3.1) 

n = I (as  in  EDUNE,  SBEACH)  or  n = 3 (as  in  CROSS)  or  other  values? 

• Different  breaking  waves  (spilling/plunging)  dissipate  their  energy  in  different 
ways,  thus  cause  different  current  patterns  and  different  sediment  transport.  How 
should  these  effects  be  represented  in  large-scale  numerical  models? 

• What  is  the  best  representation  of  the  bar-type  profile  generation  mechanism?  One 
promising  and  realistic  approach  is  to  lag  (in  the  landward  direction)  the  wave 
energy  dissipation  per  unit  volume  D, 

• What  is  the  best  representation  of  different  time  scales  of  beach  profile  evolution 
depending  on  initial  profile  and  wave  conditions? 

• What  is  the  best  representation  of  sediment  transport  under  irregular  waves? 

In  this  chapter,  the  wave  transformation  models  are  introduced  first  to  determine 
the  wave  energy  dissipation  under  both  regular  and  random  waves.  This  is  followed  by 
investigation  of  resulting  different  equilibrium  beach  profile  forms.  A new  empirical 
relationship  based  on  large-scale  experimental  data  is  then  presented  to  incorporate  the 
wave  effects  on  the  profile  scale  parameter  in  the  equilibrium  beach  profile.  Finally,  the 
detailed  development  of  the  model  and  the  numerical  solution  are  described. 


3.1  Wave  Energy  Transformation  Models 


Wave  energy  transformation  is  different  for  monochromatic  and  random  waves. 
Random  waves  typically  break  across  a much  wider  portion  of  the  beach  profile  than 
monochromatic  waves  with  the  same  average  energy,  producing  a more  evenly 
distributed  cross-shore  energy  dissipation.  Thus  they  should  be  treated  differently  in  the 

3.1.1  Wave  Breaker  Type 

Four  different  types  of  breaking  waves  have  been  identified,  namely  spilling, 
plunging,  surging  and  collapsing.  In  general,  spilling  breakors  tend  to  occur  on  beaches 
of  very  low  slope  with  waves  of  high  steepness  values;  plunging  breakers  are  associated 
with  steeper  beaches  and  waves  of  intermediate  steepness;  surging  occurs  on  high- 
gradient  beaches  with  waves  of  low  steepness;  and  collapsing  breaker  is  intermediate 
between  the  plunging  and  surging.  Galvin  (1968)  found  fairly  good  breaker-type 
correlations  with  the  well-known  Iribarren  number  (also  called  surf  similarity  parameter). 


(3.2) 


where  m is  the  beach  slope,  Ho  and  Lo  are  the  deepwater  wave  height  and  wavelength, 
respectively.  Battjes  (1974),  Okazaki  and  Sunamura  (1991),  and  Smith  and  Kraus  (1991) 


recommended  adjustments  on  the  limits  to  classify  different  breaker  types.  Table  3-1 
summarizes  die  breaker  type  classifications  from  different  resources. 


Summary  of  breaker  type  classification 


Spilling 

Plunging 

Surging 
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Different  breaker  types  decay  in  different  manners  inside  the  surf  zone.  On  steeply 
sloping  reflective  beaches,  once  waves  break  they  immediately  surge  up  the  beach  face  as 
runup,  with  little  intervening  surf.  In  contrast,  on  lower-sloping  dissipative  to 
intermediate  beaches,  there  can  be  an  extensive  surf  zone,  at  times  being  crossed  by 
several  broken  wave  bores  while  a few  small  waves  have  not  yet  broken.  Plunging 
breakers  characteristically  dissipate  their  energy  in  a concentrated  manner  in  the  region 
just  shoreward  of  the  breaker  line,  while  spilling  breakers  dissipate  their  energy  at  a 
slower,  more  uniform  rate  over  a greater  surf  zone  width.  Through  laboratory 
experiments.  Ting  and  Kirby  (1994, 1995. 1996)  found  that  turbulent  kinetic  energy  was 
transported  seaward  under  the  spilling  breaker  and  landward  under  the  plunging  breaker 
by  the  mean  flow.  Ting  and  Kirby  (1994)  stated  that 

. . .Nevertheless,  in  most  situations  there  is  a net  offshore  movement  of  sediment  in 

spilling  breakers  and  a net  onshore  movement  in  plunging  breakers. 

In  order  to  test  the  plausibility  of  their  argument,  1 4 cases  of  Large  Wave  Tank 
(LWT)  (Kraus  and  Larson,  1 988)  and  1 0 cases  of  CRIEPI  Flume  (Kajima  ct  al.,  1 982) 
experiment  data  (all  periodic  waves)  are  investigated.  Table  3-2  lists  the  test  conditions 
for  these  experiments  and  Figure  3- 1 shows  the  results.  The  criteria  of  accretive  and 
erosive  profile  based  on  the  fall  velocity  parameter  Ht/wT  and  spilling/plunging 
breaker  based  on  surf  similarity  parameter  arc  also  shown  in  the  plot. 

These  results  show  that  the  surf  similarity  parameter  alone  cannot  give  a clear  and 
correct  classification  of  profile  types,  although  the  model  is  correct  in  19  of  25  cases, 
whereas  the  fall  velocity  parameter  discriminates  correctly  in  22  of  25  cases.  The  data  in 
Figure  3-1  also  demonstrate  that  an  approximate  inverse  correlation  exists  between  the 
surf  similarity  parameter  and  the  fall  velocity  parameter  (as  discussed  by  Dean,  1987). 
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3.1.2  Wave  Transformation  Models 

3.1.2.I  Monochromatic  wave  transformation  models 

In  the  EDUNE  and  CROSS  models,  a constant  ratio  of  breaker  height  to  water 
depth  ratio  inside  the  surf  zone  (spilling  breaker  assumption)  was  used,  which  is  not 
applicable  on  non-monotonic  beach  profiles  such  as  those  containing  bar/trough  features. 
For  this  case,  energy  dissipation  occurs  only  if  dh/dy  > 0 . 

One  popular  simple  wave  model  in  the  surf  zone  is  Dally,  Dean  and  Dalrymple 
model  (Dally  ct  al.,  1984,  1985).  Hereafter,  it  will  be  referred  as  D3  model.  Other  wave 
models  such  as  those  based  on  Boussincsq  equations  are  very  time-consuming  and  not  a 
good  choice  for  very  long  lime  profile  evolution  simulation.  The  governing  equation  of 
the  D3  model  is  given  in  Eq.1.9  and  repeated  in  Eq.3.3. 


where  ECt  is  energy  flux,  ElCt  is  the  energy  flux  associated  with  the  stable  wave  that 
the  breaking  wave  is  striving  to  attain  for  the  particular  water  depth,  and  K'  is  a 
dimensionless  decay  coefficient.  For  increasing  K' , the  breaking  wave  height  dissipates 
faster  and  vice  versa.  Titus,  choosing  different  K values  in  the  DJ  model  may 
approximate  different  wave  breaking  types.  The  stable  wave  height,  H,  is  given  by 
H,  =Th  (3.4) 

where  F is  a dimensionless  coefficient.  Based  on  calibration  with  laboratory  data.  Dally 
etal.  (1985)  suggested  I'  = 0.40and  K'  =0.15  for  applications.  Dally  etal.(1985) 
presented  analytical  solutions  of  the  model  for  beaches  of  idealized  shapes  such  as 
uniform  depth,  uniform  slope  and  equilibrium  beach  profiles.  Figure  3-2  shows  an 
example  of  the  model  predicted  breaker  decay  as  compared  to  a laboratory  measurement 
for  plane  beaches  of  1/20, 1/30, 1/65  and  1/80  slopes.  The  results  show  that  the  model 
appears  to  provide  a good  representation  of  breaking  wave  decay  on  plane  beaches  of 
laboratory  scale.  The  line  H = 0.78/r'  is  also  plotted  in  the  figures  and  appears  to  be  a 
reasonable  description  of  breaker  decay  only  for  the  1/30  slope.  Figure  3-3  shows 
another  example  of  the  model  prediction  on  a prototype  scale  bar-trough  profile.  The 
simulated  wave  reaches  the  stable  criterion  in  the  deepest  portion  of  the  outermost  trough, 
shoals  on  the  inner  bar  until  the  incipient  condition  is  again  attained,  and  then  breaks 
continuously  until  the  shoreline  is  reached. 

3.1.2. 2 Random  wave  models 

There  exist  several  random  wave  models,  which  can  be  classified  into  two  groups: 
one  group  using  a representative  wave  height  with  a wave  height  distribution  assumption 
inside  the  surf  zone,  the  other  group  using  a wave-by-wave  approach. 


Figure  3-2.  Comparison  of  DJ  model  and  "0.78"  criterion  to  laboratory  data  for  I /20, 

1/30, 1/65,  and  1/80  beach  slope.  Trend  of  ( ////■')  increases  with  increasing 
(Dally  et  al..  1985).  8 P g cep  ) 
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Baltjes  and  Janssen  (1978)  first  developed  a wave  energy  dissipation  model  for 
random  waves.  The  wave  height  distribution  inside  the  surf  zone  was  assumed  to  be  a 
truncated  Rayleigh  distribution.  The  average  wave  energy  dissipation  was  modeled  as  a 
bore.  The  result  of  the  Battjcs  and  Janssen  model  is 

D = jQjpgHj 

= -f ^=-1  , H.  = 0.88*  tanhf -^-1  (3.5) 

taft  \H.)  1 0.88  J 

where  D is  die  averaged  wave  energy  dissipation.  Q,  is  the  percentage  of  breaking 
waves,  is  the  root-mean-square  wave  height,  Hm  is  a maximum  possible  wave 
height  at  water  depth  h,  f is  the  mean  wave  frequency,  * is  the  wave  number,  a and  y 
are  two  empirical  parameters.  Later,  on  the  basis  of  the  idea  of  Battjes  and  Janssen  but 
with  different  wave  height  distribution  assumptions,  Thornton  and  Guza  ( 1983),  Church 
and  Thornton  ( 1 993),  Baldock  et  al.  ( 1 998),  etc.,  have  developed  other  random  wave 
models. 

In  the  above-mentioned  models,  a narrow-banded  wave  spectrum  and  a form  of  the 
wave  height  distribution  inside  the  surf  zone  was  assumed.  However,  in  nature,  the  wave 
spectrum  is  not  always  narrow  banded  and  the  wave  height  distribution  is  very  complex, 
especially  inside  the  surf  zone.  One  method  to  circumvent  this  problom  is  the  wavc-by- 
wave  approach  (Masc  and  Iwagaki,  1982;  Dally  and  Dean,  1986;  Dally.  1992;  Larson  and 
Kraus,  1992;  Zheng,  1996;  Kuriyama  and  Nakatsukasa,  2000). 

In  the  wave-by-wave  method,  the  random  wave  field  is  represented  by  the 
collection  of  a set  of  individual  regular  waves.  The  major  assumption  inherent  in  this 
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approach  is  thal  shoaling  and  breaking  arc  noi  affcclcd  by  wave-wave  interaction.  At  any 
location,  the  probability  distribution  in  wave  height  can  be  determined  from  the  predicted 
heights  of  the  individual  waves,  and  wave  statistics  such  as  root-ntcan-square  height  are 
easily  calculated.  The  wave-by-wave  approach  has  been  verified  with  laboratory  and 
field  data  (Dally,  1992)  and  applied  to  model  the  undertow  and  longshore  current  (Larson 
and  Kraus.  1992;  Kuriyama  and  Nakatsukasa.  2000). 

Larson  (1995)  developed  an  alternative  random  wave  model  that  required 
transformation  of  only  one  representative  wave  height  without  making  any  assumptions 
about  the  pdf  of  wave  height.  Larson  (1995)  showed  that  the  model  results  were  in 
reasonable  agreement  with  the  wavc-by-wavc  approach. 

The  wavc-by-wave  approach  is  believed  to  be  more  similar  to  natural  processes  and 
is  incorporated  into  the  cross-shore  transport  model  to  simulate  the  sediment  transport 
processes  under  random  waves.  A wave  spectrum  in  the  offshore  boundary  needs  to  be 
provided  to  begin  the  calculation.  If  no  real  time  series  data  exist,  the  joint  distribution  of 
wave  period  and  heights  developed  by  Longuet-Higgins  (1983)  can  be  applied  to 
generate  wave  height  and  period  series  in  the  time  domain. 


where  R and  r are  nomialized  wave  height  and  wave  period,  respectively. 


II 


(3.7) 


m,  denotes  the  nth  moment  of  the  spectral  density,  v is  the  band-width  parameter 
determined  by  the  first  three  moments  of  the  spectrum. 


(3.8) 


The  marginal  probability  density  of  wave  height  R can  be  determined  as 

/><*)=  }p(«.r)dr  = £(i')[l  + er/(/[/v)]«exp(-/i’)  (3.10) 

where  eif(R/v)  = -^=  Jexp(-r!  )dt  is  the  error  function.  If  v = 0,  Eq.(3.10)  is  the 

Rayleigh  distribution  for  R . For  large  waves  and  small  v values,  such  as  shown  in 
Figure  3-4,  the  distribution  of  R can  be  approximated  by  the  Rayleigh  distribution. 
Therefore,  the  cumulative  probability  of  the  wave  height  becomes 

P(R)  = l-exp(-R!)  (3.11) 
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The  distribution  of  wave  period  at  fixed  values  of  the  wave  height  R is  found  by 
dividing  Eq.(3.6)  by  Eq.(3. 10).  Integrating  with  respect  to  r and  approximating  the  error 
function  erlfx)  with  tanh(1.2x),  the  cumulative  probability  of  wave  period  is  given  by 

— != i — -Li  nio.r 


i+UBh[i^ 


To  simulate  a random  wave  condition,  a random  number  Pi  uniformly  distributed 
between  0 and  I is  generated  first,  A random  wave  height,  Ri,  is  obtained  from  Eq.(3.1 1) 


R,  = yl-m~P,)  (3.13) 

Then  a second  random  number  P;  is  generated,  and  the  corresponding  wave  period,  r, , is 
calculated  from  Eq.(3.I2) 


/’.1 1 + tanh  — 


1.2«.'| 


The  dimensional  wave  height  and  wave  period  can  then  be  determined  from  Eq.  (3.7). 
Zheng  (1996)  compared  the  numerical  results  with  Longucl-Higgins's  joint  density 
distribution  for  both  narrow  ( v =0.1)  band  and  wide  band  ( v =0,1)  spectra,  and  found 
very  good  agreement  in  both  cases  (Figure  3-5). 

3.1.3  Roller  Models 

When  waves  break,  it  has  been  found  that  the  wave  energy  dissipation  does  not 
occur  immediately  at  the  initial  breaking  point  but  at  some  distance  landward  of  the 
breaker  point  (Svendscn,  1984).  It  is  within  this  distance  that  the  decay  of  organized 


wave  energy  is  convened  inlo  turbulent  kinetic  energy  that  is  not  dissipated  immediately. 
This  region  is  known  as  the  transition  region,  which  is  characterized  by  rapid  wave  decay 
and  uniform  wave  setdown,  thus  indicating  that  there  is  also  a lag  in  momentum  transfer. 
This  phenomenon  is  observed  in  both  plunging  and  spilling  breakers.  In  many 
experiments,  the  bar  crests  are  found  somewhat  shoreward  of  the  initial  breaking  points. 
Thus,  the  inclusion  of  the  transition  zone  in  a profile  model  has  an  imporiant  effect  on  the 
predicted  profiles,  particularly  near  the  bar  features  (Naim  cl  al„  1990;  Southgate  and 
Naim,  1993). 


Figure  3*5.  Comparison  of  the  numerical  results  with  the  density  function  Eq.3.37  of 

Longuot-Higgins  (Zheng,  1996).  A)  Narrow  band  (v  =0.1)  wave  spectra.  B) 
Wide  band  ( v =0.5)  wave  spectra. 

One  concept  applied  by  many  researchers  to  address  the  lag  in  energy  dissipation  is 
through  the  surface  roller  model  (Svendscn,  1984;  Dally  and  Brown,  1995;  Reniers  and 
Battjes,  1997).  The  basic  idea  of  the  roller  model  is  that  the  organized  wave  energy 
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released  at  breaking  is  first  transferred  to  the  energy  of  a roller,  which  rides  on  the  wave 
front  with  the  phase  speed  of  the  waves.  The  difference  in  velocity  between  the  roller 
and  the  water  particles  at  the  wave  front  causes  a shear  stress  between  the  roller  and 
underlying  wave  front.  The  work  done  by  the  shear  stress  equals  the  dissipation  of  roller 
energy.  The  time  needed  for  the  organized  wave  energy  to  dissipate  through  the  roller 
causes  a spatial  lag  between  the  location  of  wave  breaking  and  the  actual  dissipation 
(Naimetal..  1990). 

The  timc-avcraged  energy  balance  equation  that  contains  contributions  from  both 
the  organized  wave  motion  and  the  roller,  is  expressed  as  (Naim  et  al.,  1990;  Dally  and 
Brown,  1995) 

where  ECt  and  E,C  are  the  energy  Duxes  associated  with  the  organized  wave  motion 
and  the  roller,  respectively;  E,  and  e,  are  the  energy  and  energy  dissipation  rate  of  the 
roller,  respectively.  Svendscn  (1984)  expressed  the  roller  energy  per  unit  surface  area  in 
terms  of  the  roller  cross-section  area  A, 


The  roller  energy  dissipation  rate,  e, , is  modeled  as  the  work  performed  by  the  shear 
stress  r,  between  the  roller  and  the  organized  wave  motion. 

e,=r,C  (3.17) 

Following  Naim  et  al.  (1990),  the  shear  stress  is  deduced  by  considering  the  vertical  force 
balance  on  the  roller,  which  leads  to 


P 


(3.18) 


where/?  is  the  angle  ofthc  wave-roller  stress  vector.  Lippmann  et  al.  (1996)  showed 
that  the  model  is  insensitive  to  values  of  /?  > 10°. 

Several  formulations  for  the  roller  area  A,  have  been  proposed.  Svendsen  (1984) 
related  the  roller  area  to  the  local  wave  height  by 

A,  - OSH1  (3.19) 

Dcigaard  eta!.  (1991)  assumed  the  roller  was  similar  in  size  to  the  aerated  region  of  a 
steady,  fully  developed  hydraulic  jump  and  modeled  the  roller  area  according  to 


If  applied  in  the  transition  region,  because  of  the  dependence  on  wave  height  there  is  no 
shill  in  the  position  of  the  maximum  wave  energy  dissipation.  In  attempting  to  model  the 
shin,  Okayasu  ct  al.  (1988)  adopted  the  model 

A,=a*HL  (3.21) 

in  which  Lis  the  local  wave  length,  a and  k arc  empirically  prescribed  factors  that 
depend  upon  position  in  the  surf  zone.  The  peak  in  wave  energy  dissipation  was  shined 
landward  by  increasing  k linearly  from  zero  at  the  breakpoint  to  a maximum  value  at  the 
observed  end  of  the  transition  region  and  then  back  to  zero  at  the  shoreline.  However, 
this  method  requires  a priori  knowledge  of  the  extent  of  the  transition  region. 

Dally  and  Brown  (1995)  proposed  a di  fferent  approach  to  determine  the  roller  area. 
In  their  model,  the  rollor  area  is  not  known  in  priori,  but  obtained  by  solving  Eq.3.15 
directly  by  knowing  the  wave  heights  using  measured  data  or  from  wave  transformation 
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models.  Figure  3-6  shows  the  calculated  roller  area  using  this  method  in  comparison 
with  the  models  of  Svendsen  and  Dcigaard  ct  al.  The  shill  between  the  breakpoint  and 
the  location  where  the  roller  becomes  fully  developed  is  obvious  in  this  new  method. 
The  momentum  equation  including  roller  contribution  is 

= (3.22) 

ay  ay  ay 

where  = 2E,  = pCy  (3.23) 

This  equation  is  used  to  predict  the  wave  setup/setdown  in  the  model. 


3.2  Equilibrium  Beach  Profiles 

The  equilibrium  beach  profile  theory  developed  by  Dean  (1977)  is  based  on 
uniform  wave  energy  dissipation  per  unit  water  volume,  D. , for  a given  sediment  size 
inside  the  surf  zone. 

dF_ 
dy 


(3.24) 


(3.26) 


<3!,> 

q J.27  bock  inlo  Eq.3.24  leads  to  Eq.1.10,  which  can  be  rewritten  as 


*d.=/*P*Y 


»-y/Wik'-r,y"+5™  (3.33) 

<!j4> 


»of  r,{y) 


priori,  Ihe  developed  cross-shore  profile  model  is  applied  to  obtain  the  final  profile  and 
j’jOO.  The  equilibrium  profile  is  then  calculated  by  substituting  y,(y)  into  Eq.3.34  and 
integrating  numerically.  Figure  3-8  shows  the  initial  and  predicted  profiles,  and  Figure  3- 
9 shows  the  predicted  wave  height  to  water  depth  ratio  and  the  profile  calculated  by 
Eq.3.34.  A non-monotonic  profile  is  obtained. 


Figure  3-9.  Predicted  H/h  ratio  and  calculated  EBP  from  Eq.3.3 1 . 


3.2.4  Reevaluation  of  the  Profile  Scale  Factor 


III  tile  above  equilibrium  beach  profile  forms,  the  profile  scale  factor  A was 
assumed  to  be  only  a function  of  either  grain  diameter  or  its  equivalent  fall  velocity. 
Moore  (1982)  examined  numerous  profiles  to  quantify  the  relationship  between  A and  the 
effective  diameter  of  the  sediments  across  the  surf  zone.  Later.  Dean  (1987)  transformed 
Moore's  relationship  of  A versus  diameter  to  A versus  fall  velocity. 

A = 0.0671V1’"  (3.35) 

Table  3-3  lists  the  recommended  A values  versus  the  median  sediment  size  (reproduced 
from  Dean  and  Dalrymple,  2002). 


Thus,  this  representation  cannot  explain  one  phenomenon  that  occurs  in  natural  beaches: 
the  profiles  are  of  milder  slope  in  the  winter,  compared  to  the  summer  conditions 
associated  with  milder  and  less  steep  waves.  The  fall  velocity  parameter  includes 
information  of  both  sediment  size  and  wave  conditions.  The  hypothesis  tested  here  is 
that  there  exists  a relationship  between  the  profile  scale  factor  A and  the  fall  velocity 


(3.36) 


67 

where  the  first  term  is  determined  from  Table  3-3,  and  the  second  term  is  a modification 
term  which  is  a function  of  the  fall  velocity  parameter.  To  investigate  this  modification 
function,  the  large-scale  wave  flume  experimental  data  from  LWT  (Kraus  and  Larson, 
1988),  CRIEPI  (Larson  and  Kraus,  1989),LWF,96/97(Detteetal„  1998)  and  LIP1ID 
(Roelvink  and  Rcnicrs,  1995)  are  analyzed.  The  detailed  descriptions  of  these 
experiments  are  included  in  the  following  chapters  for  model  calibration  and  verification 
purposes.  Among  them,  LWT  and  CRIEPI  used  monochromatic  waves,  while  the 
experiments  in  LWF’96/97  and  LIP1 1 D were  conducted  with  random  waves. 

The  analysis  procedure  is  as  follows:  first  using  the  DJ  model  or  from  the 
measurements,  determine  the  breakpoints  and  the  corresponding  breaking  wave  heights 
(significant  wave  heights  are  used  in  the  case  of  random  wave  conditions)  in  the  final 
measured  profiles;  then,  applying  the  least-square  method  from  the  shoreline  positions  up 
to  the  breakpoint  positions,  determine  the  best-fit  A values  by  assuming  that  the  final 
profiles  are  in  equilibrium  (although  some  of  them  may  not  be).  Figure  3-10  shows  one 
example  from  CE300  in  the  LWT  experiments.  A plot  of  the  fall  velocity  parameter  and 
the  ratio  of  the  best-fit  A values  to  the  recommended  values  in  Table  3-3  are  presented  in 

Upon  further  investigation  of  the  experimental  data,  the  following  arc  believed  to 
be  partially  responsible  for  a major  portion  of  the  scatter:  ( I ) there  may  be  an  effect  the 
initial  slope  on  the  final  equilibrium  profiles,  the  initial  slope  ranges  from  1:50  to  1:10; 
(2)  some  profiles  have  not  reached  equilibrium,  particularly  those  lower  data  points  on 
the  left  in  Figure  3-11,  because  these  cases  all  have  initial  mild  slopes  under  small  wave 


action  and  the  test  durations  are  short  (the  time  scale  of  profile  accretion  is  much  longer 
than  erosion);  (3)  the  final  profiles  have  multiple  bars  and  a wide  trough  landward  of  the 
primary  bar  in  tests  such  as  CE500,  CE501 , CRIEP13-4,  CRIEPI4-3  and  CR1EPI6- 1 . An 
attempt  is  made  to  fit  these  profiles  just  to  the  second  breakpoint  near  the  inner  bats,  with 
better  results  than  the  original  (Figure  3-12).  Figure  3-1 3 shows  an  example  of  the 
adjustment  in  CESOO. 


3 

I 


Figure  3-10.  Best-fit  equilibrium  beach  profile  for  CE300,  The  A value  based  on  Table  3- 
2 is  0.106m  . 


Figure  3-13.  Modified  best-fit  equilibrium  profile  for  CE500.  The  A value  based  on 
Table  3-2  is  0.106  mf'3. 

Figure  3-14  shows  Ihe  relationship  based  only  on  the  LWT  profiles  (all  commenced 
with  an  initial  slope  of  1 : 1 5)  and  most  of  the  final  profiles  in  this  experiment  had  almost 
achieved  equilibrium.  As  Figure  3-14  illustrates,  a reasonable  fit  to  the  data  is 


(3.37) 
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This  equation  will  be  included  in  the  mode!  to  modify  the  profile  factor.  A,  in  the 
determination  of  uniform  wave  energy  dissipation  per  unit  volume,  D. . Of  special 
interest,  seasonal  profile  changes  can  be  interpreted  using  Eq.3.37:  from  winter  to 
summer,  the  A value  becomes  larger  with  the  fall  velocity  parameter  decreasing,  thus  the 
profile  becomes  steeper  in  summer  and  the  shoreline  advances  seaward. 


3.3  Cross-Shore  Profile  Evolution  Model 


As  discussed  in  Dean  and  Dalrymple  (2002),  an  equilibrium  profile  represents  a 
balance  of  destructive  and  constructive  forces  noting  on  the  beach.  Of  those  forces 
recognized,  the  dominant  destructive  force  inside  the  surf  zone  is  the  high  turbulence 
level  caused  by  wave  breaking.  The  transport  models  based  on  wave  energy  dissipation 
and  equilibrium  beach  profile  concepts,  such  as  EDUNE,  SBEACH  and  CROSS,  capture 
this  force.  Although  these  models  lack  detailed  descriptions  of  the  physical  processes 
related  to  sediment  transport,  they  have  yielded  reasonable,  useful  and  encouraging 
results,  particularly  for  storm  wave  conditions  which  cause  beach  and  dune  erosion. 
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Thus,  ihc  model  developed  in  this  research  uses  the  same  strategy  with  some 
modifications  made  based  on  recent  and  additional  analysis. 

33.1  Global  Consideration  of  Transport  Rate  Exponents 

In  the  transport  rate  relationship  (Eq.3.1),  two  different  exponent  (n)  values  have 
been  used:  n =1  (EDUNE  and  SBEACH),  and  n=3  (CROSS).  The  experimental  data 
from  LWT  (Kraus  and  Larson,  1988)  arc  examined  to  investigate  this  relationship. 

In  the  case  of  n = 1 , Eq.3. 1 can  be  integrated  from  the  shoreline  position  y„  (where 
the  wave  energy  is  assumed  to  be  close  to  zero)  to  the  wave  breaking  point  yb  (assuming 
no  wave  energy  dissipation  outside  the  surf  zone), 

\hQdy  = K j (AD  - AD.  )dy  (3.38) 

This  leads  to  Eq.3.39,  die  relationship  between  the  total  wave  energy  entering  the  surf 
zone  and  the  total  sediment  transport  (determined  from  the  profile  measurements). 


Figure  3-15  shows  the  results  for  case  CE300,  CE400  and  CE40I,  where  the 
calculation  has  been  carried  out  at  the  times  the  measured  profiles  are  available. 
However,  for  n = 3,  direct  integration  is  impossible  unless  the  distribution  of  wave 
energy  dissipation  is  known.  Thus,  the  D3  model  is  applied  to  predict  the  wave  energy 
transfomiauon  within  Ihc  surf  zone.  Eq.3.1  is  then  integrated  numerically.  The  results 
are  presented  in  Figures  3-16  and  3-17  for  n - 1 and  n = 3,  respectively. 
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From  the  above  comparisons,  n = I appears  to  be  the  better  representation  of  the 
relationship  between  transport  rates  and  wave  energy  dissipation.  However,  from 
physical  considerations,  n - 3 is  considered  to  be  consistent  with  both  Froude  number 
and  the  fall  velocity  parameter  modeling  similarity  criteria  (Zheng,  1996;  Zheng  and 
Dean,  1 997).  Some  discrepancies  may  be  due  to  the  inaccuracy  of  modeling  wave  energy 
dissipation  and  determination  of  transport  rates  from  measured  profiles.  Thus,  both  of 
these  two  n values  are  included  in  the  model  for  further  investigation. 

3.3.2  Transport  Rates 

A modified  Eq.3.1  that  includes  the  gravity  effect  is  used  to  obtain  the  transport 
rates  inside  the  surf  zone. 

o.<i) 

where  a is  a gravity  coefficient,  which  equals  to  D.  I m according  to  Dean  (1991),  m is 
the  beach  face  slope.  Three  different  wave  energy  dissipation  models  have  been 
examined  and  compared  in  Section  3.2;  spilling  breaker  assumption,  D model  and  roller 
model.  The  first  two  models  cannot  treat  the  transition  zone  directly  and  the  resulted 
equilibrium  profiles  are  monotonic.  In  the  program  SBEACH,  the  transport  rates  in  the 
transition  zone  are  considered  to  decay  exponentially  from  the  plunge  point  to  the  break 
point.  In  contrast,  Lee  et  al.  (1996)  assumed  a Gaussian  distribution  in  that  region.  On 
the  other  hand,  the  roller  model,  especially  the  formula  ofDally  and  Brown  (1995),  can 
directly  account  for  the  lag  of  energy  dissipation  and  a non-monotonic  equilibrium  profile 
can  be  reached.  Thus,  theD3  model  coupled  with  the  roller  model  ofDally  and  Brown  is 
chosen  in  the  present  study  to  predict  the  wave  height  transformation  and  energy 
dissipation  inside  the  surf  zone. 
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333  Numerical  Modeling  Process 
333.1  Grid  system 

Two  types  of  grid  representations  of  the  physical  domain  have  been  used  for  cross- 
shore  numerical  modeling  as  shown  in  Figure  3-18.  In  the  first  type  shown  in  Figure  3- 
1 8A,  the  cells  are  finite  increments  of  the  distance  variable,  y.  whereas  in  the  second  type 
shown  in  Figure  3-I8B,  the  computational  cells  are  formed  by  finite  increments  of  the 


depth  variable,  h.  The  advantage  of  the  first  type  is  that  it  can  easily  represent  bar 
profiles,  thus  it  is  chosen  in  the  present  mode]  with  a staggered  representation  such  that 


the  depth,  h,  is  located  at  the  center  of  each  cell  and  the  transport,  Q,  is  at  the  grid 
boundaries.  H and  C,  are  at  the  Q-points,  while  D is  at  the  h-points.  The  cross-shore 


coordinate,  y,  is  directed  offshore,  thus  positive  Q indicates  offshore  transport. 


3.3.3. 2 Governing  equations 


The  governing  equations  in  the  model  are  listed  below  for  convenient  reference: 
• D'  Model 


SEC,  K r , 

— = 7— ^[CC,  - EC  \ 
By  h + tj 1 ' 

■ Dally  and  Brown's  roller  model 


• Momentum  equation 


S„,„  = S„.,  =2 E,  = pC  — 


(3.42) 


(3.43) 


(3.44) 


(3.45) 


(3.46) 


where  h includes  wave  setup/setdown,  which  is  calculated  by  iteration. 

3.3.3.3  Numerical  solution  scheme 

There  are  two  types  of  numerical  solutions  available:  the  “explicit"  and  "implicit" 
methods.  Tile  explicit  method  solves  the  transport  and  continuity  equations  sequentially 
whereas  the  implicit  method  solves  these  two  equations  simultaneously.  The  primary 


advantage  of  the  explicit  method  is  simplicity  in  formulation  and  programming,  and  also 
the  ease  of  incoiporating  different  transport  rate  formulations;  however,  in  order  to  avoid 
numerical  instabilities,  very  small  time  steps  can  be  required.  The  implicit  method 
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allows  much  larger  time  steps  while  maintaining  the  solution  stable.  The  implicit  scheme 
is  employed  to  solve  the  governing  equations  in  the  present  model.  The  finite  difference 
equations  for  each  module  arc  presented  below.  In  these  equations,  the  index  i is  used  to 
denote  the  number  of  a specific  grid  point  and  increases  offshore.  Linear  wave  theory 
and  shallow  water  assumptions  are  applied  to  calculate  wave  related  quantities. 

For  the  wave  module,  the  numerical  calculations  start  at  the  seaward  end  of  the  grid 
and  proceed  onshore.  Wave  setup  is  included  in  the  water  depth  as  an  unknown,  thus  an 
iteration  procedure  is  needed  until  the  wave  setup  calculation  meets  certain  accuracy 
criteria.  From  a specific  Q-poinl  to  the  next,  Eqs.3.42  to  3.44  are  solved  numerically  one 
after  another.  The  difference  equation  for  Eq.3.42  is 

(£C,)W -(EC,),  A'Tl  1 

Ay  °7^1lt(£C,>'"  +(£C')'1_(£‘C»)'J  (3‘47) 


whcre  (£Cp),  =^yg[rh,f C,i  *C*'"  (3.48) 

and  Ay  is  the  space  step.  From  Eq.3.47,  the  energy  flux  at  point  i is, 

(EC,  ),  = [(£C„ )„,(!  -B,)  + B,(EC„  ),  ] (3.49) 


After  the  energy  (lux  has  been  calculated,  the  corresponding  wave  height  is  determined  as 

""bid 


The  difference  equation  for  the  roller  model  (Eq.3.43)  is  calculated  to  determine  the 
roller  area. 


-(■4,C1),]=6IK1.1  H 


Using  the  wave  height  and  the  roller  area,  wave  setup/setdown  is  given  from  Eq.3,44, 
expressed  in  finite  difference  form  as 


Once  the  desired  accuracy  is  reached,  the  energy  dissipation  per  unit  volume  for  the  roller 
model  at  each  Q-point  is  determined  by 


The  transport  rate  distribution  within  the  surf  zone  is  calculated  from  Eq.3.45, 
written  in  difference  form  at  the(m+l)th  time  step, 


er*1=*[4""-D.+t 

Rewriting 


The  roller  area  is  assumed  to  be  constant  over  a time  step.  From  Eq.3.55, 


(3.57) 

(3.58) 


Thus,  Eq.3.56  becomes. 


a"' + -a  *sH 


a; =1.0 

^-4— “^"(S*i) 
*-<* 


a‘*-4w'-a“£?H£-5)""‘(S*sH 

(3.62) 


whliea~*-<r±er. 


MT'  + BfilC+CjZZ-G,  (3.64) 

«.='•«.  c,=-a„  a, =^(a:,  -sr) 


» + l£i: 


(3.66) 


3. 3.3. 4  Wave  runup 


The  uprush  limit  is  based  on  the  wave  runup.  Hunt  (1958)  proposed  an  equation 
for  the  wave  runup  as  a function  of  bed  slope,  wave  height  and  wave  steepness  in  the 
study  of  seawall  and  breakwater  design,  expressed  as 


where  R is  the  wave  runup  height,  tan  p is  the  average  bed  slope  from  the  runup  limit  to 
the  incipient  breaking  point. 

3.3.3.5  Limiting  onshore  and  offshore  slopes 

The  program  adopts  die  idea  of  CROSS  (Zheng,  1 996;  Zheng  and  Dean,  1 997)  and 
allows  specification  of  maximum  onshore  and  offshore  slopes.  These  limitations  are 
required  at  the  upper  limit  of  runup  where  seaward  or  landward  transport  can  cause  the 
angle  of  repose  of  the  profile  to  be  exceeded  above  this  elevation.  Slopes  exceeding  the 
allowed  can  also  occur  at  the  seaward  end  of  the  computational  domain  where  the  most 
seaward  active  contour  is  displaced  seaward,  thus  creating  a locally  sleep  slope.  At  each 
time  step,  the  program  checks  and  if  it  is  determined  that  a local  slope  exceeds  the 
specified  values,  avalanching  occurs  with  the  sand  redistributed  to  the  adjacent  contours 
until  the  maximum  allowable  slope  is  reached.  The  allowable  slopes  from  Zheng  (1996) 
are  taken  as  0.18  and  0.22  for  onshore  and  offshore  slope,  respectively,  based  on  large 
wave  flume  data. 

3. 3.3. 6 Random  waves 

hi  the  case  of  random  wave  conditions,  the  wave-by-wave  approach  is  applied  in 
the  model.  The  wave  heights  and  corresponding  wave  periods  are  generated  as  described 


4'AC + B]vr +c>;-  = g;  (3.69) 
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CHAPTER  4 

MODEL  CALIBRATION  WITH  LABORATORY  EXPERIMENTS 


4.1  Introduction 


Application  of  the  cross-shore  model  discussed  in  Chapter  3 requires  quantification 
of  several  coefficients.  To  reduce  the  number  of  free  variables  and  facilitate  the 
calibration  process,  the  coefficients  in  the  wave  and  roller  model,  namely  decay 
coefficient  K , stable  wave  coefficient  r , and  roller  dissipation  coefficient  pd , are  taken 
as  the  suggested  values  at  this  stage,  that  is,  AT'  = 0.15,  T = 0.4 . and  pa  =0.1.  Thus, 
two  coefficients  are  to  be  determined,  the  exponent,  n,  and  the  transport  coefficient,  K. 
The  linear  (n  = I ) and  nonlinear  (n  = 2, 3)  transport  relationships  are  used  for  comparison 
of  model  prediction  with  several  large-scale  laboratory  experiments.  The  similarities  for 
these  models  are  that  they  should  converge  to  the  same  target  profile  for  unchanging 
wave  conditions  and  unlimited  time.  The  only  difference  is  that  the  evolution  with  larger 
n value  will  evolve  more  rapidly  at  the  beginning  and  more  slowly  later  than  the 
relationship  with  smaller  n values. 

As  an  objective  criterion  forjudging  agreement  between  the  simulated  and 
measured  profile,  the  root-mean-square  of  the  difference  between  measured  and 
calculated  depths  is  used. 


in  which  superscripts  m and  p refer  to  measured  and  predicted  profiles,  respectively,  and 
N is  the  number  of  data  points.  The  errots  between  the  predicted  and  observed  profile 


depths  are  obtained  at  the  times  for  which  measured  profile  information  is  available.  The 


best-fit  K value  is  determined  as  the  value  yielding  the  overall  least  RMS  error.  In  order 
to  avoid  scale  effects  on  profile  evolution,  only  large-scale  wave  flume  experiments  are 
used.  Dctte  et  al.  (2002)  provided  a detailed  overview  of  the  existing  large  wave  flumes 
in  the  world  and  the  profile  evolution  related  experiments  available  from  these  flumes. 
Three  different  laboratory  data  sets  were  available  and  are  examined  in  this  study  for  the 
present  model  calibration:  "Large  Wave  Tank”  (LWT)  experiments  in  1956-1957  and 
1962  by  Saville  (data  are  from  Kraus  and  Larson.  1988);  the  Large  Wave  Flume  of  the 
Central  Research  Institute  of  Electric  Power  Industry  (CRIEPI)  in  1982  (Kaj  ima  et  al. 
1982;  Larson  and  Kraus.  1989);  and  German  wave  flume  experiments  (GWK)  in 
Hannover,  Germany,  in  1986/1987  (Detle  and  Uliezka,  1987)  and  inl996/97  as  part  of 
the  MAST  111-SAFE  project  (Dette  et  al..  1998).  Among  them,  LWT  and  CRIEPI  were 
carried  out  with  monochromatic  waves;  GWK  had  both  monochromatic  and  random 

4.2  LWT  Experiments 
4.2.1  Description  of  Experiments 

Two  series  of  movable-bed  model  experiments  on  beach  profile  change  were 
conducted  in  the  lank  of  the  Beach  Erosion  Board  1956/57  and  1962  (reported  by  Kraus 
and  Larson,  1988).  Tile  two  series  differ  mainly  in  the  grain  size  of  the  bed  material 
used.  Nine  major  tests  were  run  with  each  of  two  sand  sizes  using  essentially  the  same 
water  depth  and  wave  conditions.  Tests  for  a given  grain  size  were  characterized  by 
waves  of  varying  steepness,  producing  either  accreting  or  eroding  profiles.  All  tests  were 
with  monochromatic  waves.  The  experiments  included  extensive  profile  surveys  during 
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each  test  to  document  the  approach  to  equilibrium.  Only  limited  hydraulic  measurements 
were  made,  mainly  wave  heights  at  incipient  breaking  and  run-up  heights. 

The  experiments  were  conducted  when  the  lank  was  located  at  Dalccarlia 
Reservation  in  Washington.  DC.  The  concrete  tank  was  194  m long.  4.6  m wide,  and  6.1 
m deep.  The  initial  profile  slope  for  all  tests  but  two  was  1 : 1 5.  The  wave  periods  ranged 
between  3.75  and  1 6.0  s and  generated  wave  heights  between  0.55  and  1 .68  m in  the 
horizontal  part  of  the  tank.  The  wave  height  and  period  were  held  constant  except  in  one 
case  when  the  water  depth  was  varied  in  order  to  simulate  a tide  (Test  911).  The  water 
depths  ranged  from  3.5  to  4.6  m and  two  grain  sizes  were  used  with  median  diameters  of 
0.22  mm  in  1956/57  and  0.4  mm  in  1962,  respectively.  Tests  510  and  610  were 
conducted  fiom  an  initial  profile  that  was  non-planar  and  Test  1 10  was  a repetition  run  of 
Test  100.  Of  the  15  generic  cases  studied  (plane  initial  slope,  different  wave  and 
sediment  conditions),  eight  terminated  with  crosional  profiles  (Tests  100, 300, 400, 500, 
700, 401, 501,  and  901)  and  seven  with  accretionary  profiles  (Tests  200, 600, 101, 201, 
301, 701,  and  801). 

In  Test  100,  the  profile  eroded  back  to  the  end  of  the  tank.  The  water  level  in  Test 
700  varied  during  the  experiment.  In  all  the  accretionary  cases,  a majority  of  sand  was 
carried  onshore  from  offshore  of  the  surf  zone,  however,  the  model  calculation  domain  is 
ended  at  the  breaking  point,  thus  cannot  predict  these  onshore  transports.  Thus,  only  five 
crosional  cases  (Tests  300, 400, 500, 401  and  501 ) are  used  for  the  model  calibration. 
Table  4- 1 summarizes  the  conditions  of  the  tests  and  the  corresponding  fall  velocity 
parameters  based  on  the  breaking  wave  heights. 


4.2.2  Modeling  Results 


The  best-fit  K values  for  both  linear  (n  — 1 ) and  nonlinear  (n  =2, 3)  transport 
relationships  of  the  five  LWT  tests  arc  listed  in  Table  4-2.  It  appears  that  transport 
relationship  with  n = 2 provides  overall  least  rool-mcan-square  errors  and  best  fits  to  the 
observed  profile  evolution. 

The  comparisons  of  predicted  and  measured  profile  evolutions  are  presented  in 
Figure  4-1  to  Figure  4-5  for  the  five  cases,  respectively.  In  each  case,  the  profiles  at  four 
different  times  are  compared.  A primary  breaking  point  bar  is  predicted  at  each  case,  but 
no  secondary  inshore  bar  is  generated.  Also  the  model  cannot  predict  the  wide  trough 
area  after  the  primary  olTshore  bars.  The  nonlinear  relationships  predict  a faster  profile 
evolution  at  the  beginning  and  slower  later,  as  the  measured  profiles  show  the  same  trend. 


Test 

No. 

Grain  size 

Wave 
height  (ml 

Wave 
period  (s) 

depth  (m) 

Duration 

(hrs) 

Ht 

300 

400 

500 

0.22 

1.68 

11.33 

4.27 

50 

5.53 

1.62 

5.6 

4.42 

40 

10.0 

1.52 

3,75 

4.57 

100 

12.27 

SOI 

0.40 

1.62 

5.6 

4.42 

66 

536 

1.52 

3.75 

4.57 

60 

6.57 

fable  4-2.  Best-fit  K values  for  LWT  experiments. 

Test 

No. 

n=l 

n = 2 

n = 3 

Best-fit  K 
(m4/N) 

Best-fit  K 
(m6s/N2) 

Best-fit  K 
(m8s2/NJ) 

300 

13x10* 

0.361 

0.79x10* 

0.58x1 0*'° 

0.348 

400 

2.0x1  O'6 

0.250 

1.4x10* 

0.203 

0.99x10*'° 

0.222 

500 

1.0x10* 

0.381 

1.1x10* 

0.365 

1.1x10'° 

0.393 

401 

0.83x10* 

0.267 

0.47x10* 

0.261 

0.27x10*'° 

0.270 

501 

0.98x10* 

0.303 

0.71x10* 

0.306 

0.62x10*'" 

0.319 

4.3  CRIEPI  Wave  Flume  Experiments 

The  wave  flume  is  operated  by  the  Central  Research  Institute  of  Electric  Power 
Industry  in  Japan  between  1979  and  1983  (Kajima  et  al.,  1982),  It  is  made  of  reinforced 
concrete,  and  is  205  m long,  3.4  m wide,  and  6 m deep  in  the  1 15  m-Iong  upwave  flat 
section,  which  is  joined  to  a shoreward  section  of  1/15  slope.  As  in  the  LWT 
experiments,  two  different  median  grain  sizes  were  employed,  although  slightly  coarser 
in  the  CRIEPI  experiments  (0.27  and  0.47  mm).  A majority  of  the  tests  started  from  a 
plane  beach  slope,  but  the  initial  slope  was  varied,  ranging  from  1 :50  to  1:10.  Both 
erosional  and  accretionary  cases  were  investigated.  Only  monochromatic  waves  were 
used  in  the  experiments  keeping  the  water  level  fixed  during  a test.  Plunging,  spilling, 
and  surging  breaking  waves  were  observed,  although  plunging  breakers  occurred  in  the 
majority  of  cases.  The  experimental  program  consisted  of  24  monochromatic  tests  with 
wave  periods  ranging  between  3.0  and  12.0  s and  generated  wave  heights  between  0.3 
and  1.8  m. 

Extensive  surveys  were  carried  out  during  the  experiments  at  a very  high  spatial 
resolution  (al  least  every  50  mm,  sometimes  every  mm),  whioh  allowed  for  the  study  of 
micro-scale  features  such  as  ripples.  Larson  and  Kraus  (1989)  used  the  CRIEPI  data  to 
quantify  and  model  beach  profile  evolution.  Of  the  1 7 tests  identified  as  giving  distinct 
responses  corresponding  to  bar  or  berm  profiles.  Tests  1-8, 2-1, 3-1, 3-2, 3-4, 4-2, 4-3, 5- 
2, 6-1.  and  6-2  produced  bar  profiles  and  Tests  1-1, 1-3, 2-2, 2-3, 3-3, 4-1,  and  5-1 
produced  berm  profiles.  Here,  Test  6-2  with  a 0.27  mm  median  grain  size  is  selected  to 
calibrate  the  model  and  the  results  are  presented  in  Table  4-3.  Figure  4-6  shows  the 
comparisons  between  the  measured  and  the  best-fit  predicted  profiles  for  four  time  scries. 
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Table  4-3.  Bcsl-fil  K values  for  CRIEPI  experiment. 


No. 

Test  cond. 

n-1 

n = 2 

n = 3 

H 

T 

Best-fit  K 
(mJ/N) 

Best-fit  K 
(m0s/N;) 

Best-fit  K 
(mV/NJ) 

6-2 

1.12 

7.5 

2.3x10'*’ 

0.186 

1.2x10'® 

0.181 

0.7xl010 

0.201 

4.4  German  Large  Wave  Flume 

The  German  "large  wave  flume"  in  Hannover  is  324  m long,  7 m deep  and  5 m 
wide.  Since  1986,  a series  of  beach  profile  experiments  have  been  conducted  in  the 
flume.  Only  tests  from  1986/1987  and  1996/1997  were  available  for  the  present  study. 
4.4.1  1986/1987  Experiments 

The  purpose  of  the  first  set  of  experiments  in  1 986/1987  was  to  study  beach  and 
dune  erosion  under  monochromatic  and  random  waves  (Dettc  and  Uliczka,  1987).  Two 
of  these  experiments  were  conducted  with  the  same  constant  wave  conditions  and 
different  initial  profiles.  Regular  waves  with  a wave  height  of  1.5  m and  a period  of  6 s 
were  generated  in  a water  depth  of  5 m.  The  sand  used  for  both  experiments  had  a mean 
diameter  of  0.33  mm,  which  corresponds  to  a fall  velocity  of  5 cm/s.  Two  initial  profiles, 
"dune  without  foreshore"  and  "dune  with  foreshore",  were  used.  The  "dune  without 
foreshore"  had  a dune  crest  of  2 m above  still  water  level  and  a seaward  slope  of  1 :4 
down  to  the  channel  floor.  The  "dune  with  foreshore"  had  a slope  of  1 :4  from  the  dune 
crest  of  2 m above  still  water  level  to  I m water  depth  followed  by  a slope  of  1 :20  down 
to  the  channel  floor. 

Here,  the  case  "dune  without  foreshore"  is  used  for  model  calibration.  Table  4-4 
provides  the  results,  and  Figure  4-7  presents  the  comparisons  between  predicted  and 
measured  profiles  at  four  different  time  series. 
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Test 

n = 1 

n = 2 

n = 3 

Best-fit  K 
(m*/N) 

Best-fit  K. 
(m6s/N:) 

Best-fit  K 
(mV/N!) 

Dune  w/o 
foreshore 

2.3X10-* 

0.379 

0.79x10'* 

0.313 

1.5x10'" 

4.4.2  1996/1997  Experiments 

From  November  lo  December  1 996,  February  lo  March  1 997  and  June  lo  Augusl 
1997.  in  conjunciion  wilh  Ihc  MAST  III-SAFE  project,  experiments  were  carried  out  in 
the  Grosser  Wellenkanal  in  Hannover  to  study  the  development  of  the  underwater  profile 
and  sand  losses  from  the  beach  by  various  water  level  and  wave  conditions  (Dette  et  al., 
1998;  Dette  et  at.,  2002).  The  tests  were  divided  into  two  major  phases  focusing  on 
beach  and  dune  stability  under  given  conditions  wilh  variable  sea  levels. 

In  Test  Phase  I,  attention  was  concentrated  on  beach  stability  without  structural  aid.  The 
beach  was  formed  from  sands  with  a median  diameter  of  0.3  mm  and  a fall  velocity  of 
0.042  m/s.  The  underwater  profile  from  4 m above  the  flume  floor  was  shaped  to  the 
equmonum  profile  h = 0.12y"\  The  beach  above  the  still  water  level  (4  m above  the 
flume  floor)  had  slopes  of  1:20. 1:15, 1:10,  and  1:5.  The  initial  profile  geometries  of  Test 
Phase  I are  shown  in  Figure  4-8.  Each  test  slope  was  first  subjected  to  normal  wave  and 
water  level  conditions  using  a TMA  spectrum  with  HM  = 0.65  m and  7.  = 5 s.  This 
initial  test  lasted  until  the  profile  reached  an  equilibrium  state  for  the  selected  wave 
parameters,  and  reused  as  the  initial  profile  for  the  profile  change  due  to  storm  surge 
conditions,  by  raising  the  water  level  by  1 m (channel  floor  +5  m)  and  generating  storm 
waves  (TMA  spectrum  with  Hm0  =1.2  m and  T„  = 5 s).  The  peak  period  for  the 
generated  TMA  spectra  was  6.6  s (normal  and  storm  conditions).  The  spectral  length  was 
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1 5 min  and  consisted  of  1 65  waves.  Figure  4-9  presents  the  initial  and  final  measured 
profiles  under  storm  wave  condition  in  Test  Phase  I. 

The  model  calibration  uses  the  wave-by-wavc  random  wave  model  discussed  in 
Chapter  3.  The  best-fit  K values  for  die  four  tests  under  storm  wave  conditions  are 
summarized  in  Table  4-5  and  the  comparisons  between  predicted  and  measured  profile 
evolutions  are  shown  in  Figure  4-10  to  Figure  4-13,  respectively.  The  developed  bars 
near  the  normal  water  line  (4  m above  the  floor)  are  not  generated  by  the  model.  The 


1996/1997  experiments. 


Test  No. 

n = 1 

n = 2 

n-3 

Best-fit  K 
(m4/N) 

(m'VN;) 

Best-fit  K 

A9  (1:20) 

0.99x10'" 

0.0676 

0.68x1  O'8 

0.0728 

0.47x1  O'10 

0.0758 

B2  (1:10) 

2.5x1  O'" 

0.0765 

l.OxIO"8 

0.0828 

I.5xl0'° 

0.0883 

C2  (1:5) 

I.7X10" 

0.0721 

O.SlxlO'8 

0.0834 

0.75x1  O'10 

0.0919 

H2  (1:15) 

UxlO-" 

0.0662 

0.66x1  O'8 

0.0765 

0.36x10''° 

0.0841 
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4.5  Summary  of  Laboratory  Calibration 
A total  of  eleven  large-scale  experiments  from  three  different  wave  flumes  are 
employed  for  the  model  calibration  and  comparison  of  the  linear  and  nonlinear  transport 
relationships.  Seven  of  these  tests  were  conducted  with  monochromatic  wave  conditions 
and  the  other  four  were  carried  out  with  random  waves.  These  tests  arc  characterized 
with  different  sediment  sizes,  initial  profiles,  and  incident  wave  conditions. 

The  minimum  root-mean-squarc  errors  between  predicted  and  measured  profiles 
arc  used  as  the  objective  in  the  calibration  procedure.  The  experimental  conditions  and 
the  best-fit  results  are  summarized  in  Table  4-6.  Note  that  the  best-fit  K values  obtained 
are  based  on  total  time  series  available  for  each  individual  test , not  just  one  particular 
measured  profile.  The  average  best-fit  K value  of  the  linear  relationship  (n  = 1 ) is 
1 .72x1  O'6  m4/N,  with  a variation  from  -52%  to  133%,  whereas  the  best-fit  K values  range 
from  -49%  to  72%  of  the  average  value  of  0.93x  1 O'®  m°s/N2  for  n = 2 and  from  -78%  to 
121%  of  the  average  value  of  0.68x10''°  m8s2/N!  for  n = 3,  respectively.  The  nonlinear 
relationship  with  n - 2 has  a narrower  range  than  that  for  n = 1 and  3.  Of  the  eleven  tests, 
the  nonlinear  relationship  with  n = 2 has  the  best  fit  for  six  of  seven  monochromatic  wave 


cases,  while  n = 1 provides  the  best  fit  for  four  random  wave  cases. 

To  further  investigate  the  transport  relationship,  the  best-fit  K values  (shown  in 
Figures  4-14  to  4-1 7)  for  each  individual  measured  profile  are  determined  using  the  same 
criteria  for  LWT  tests  300, 400, 500  and  CRIEPI  lest  6-2.  In  most  cases,  the  best-fit  K 
values  for  n = I decrease  with  lime,  and  increase  with  time  for  n = 3,  but  have  no 
apparent  trend  for  n = 2.  If  the  transport  coefficient  in  the  model  is  assumed  to  be  a 


constant,  the  model  with  n - 2 can  have  a better  simulation  of  tile  profile  change  history. 
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Figure  4-15.  Besl-fit  K values  with  time  for  LWT  Test  400  with  n = 1,2  and  3. 
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Another  parameter  characterizing  model  skill  is  the  profile  foreshore  volumetric 
erosion  (Zheng,  1 996,  Zheng  and  Dean,  1 997).  Figure  4- 1 8 shows  comparisons  of  the 
measured  and  the  best-fit  predicted  eroded  volume  for  the  above-investigated  tests.  The 
best-fit  RMS  errors  are  presented  in  Table  4-7.  The  transport  relationship  with  n = 2 
gives  the  best  fit  for  Tests  CE300  and  CE400,  while  n = 1 has  the  best  lit  for  Tests 
CE500  and  CRIEPI6-2.  One  unusual  feature  in  CE500  case  is  that  there  is  a rapid  change 
in  the  eroded  volume  after  60  houts  wave  action.  If  excluding  the  profile  evolution  after 
60  hours,  n = 2 gives  the  best-fit  eroded  volume  for  this  lest. 

Based  on  these  results,  it  is  determined  that  the  nonlinear  transport  relationship 
with  n = 2 is  more  appropriate  in  the  cross-shore  profile  model.  The  transport  coefficient 
K = 0.93'  1 O'"  m"s/N  is  adopted  according  to  the  average  of  the  best-fit  K values  in  the 
eleven  tests. 


Table  4-7.  Best-fit  RMS  errorsof  eroded  volume  (unit:  nr). 


• The  measurements  after  60  hours  arc  excluded  in  the  calculation, 

**  The  RMS  error  between  predicted  and  measured  eroded  volume,  defined  as  Eq.4. 1 
except  change  profile  depth  to  eroded  volume. 


CHAPTERS 

MODEL  EVALUATION  AND  SENSITIVITY  ANALYSIS 
In  Chapter  4 eleven  large-scale  wave  flume  laboratory  tests  were  selected  to 
calibrate  the  developed  cross-shore  profile  model.  Based  on  the  calibration  results,  the 
nonlinear  transport  relationship  with  n = 2 was  found  to  yield  the  overall  least  root-mean- 
squarc  error  for  profile  evolution.  The  average  best- fit  K value  of  0.93x1  O'8  mVN1  was 
recommended  for  model  application.  The  purpose  of  this  chapter  is  to  evaluate  the  cross- 
shore model  with  the  recommended  average  K value.  First,  the  predicted  profile 
evolutions  with  the  average  K value  are  compared  with  the  results  of  the  individual  best- 
fit  K values  using  the  same  laboratory  data  in  Chapter  4.  Then  a sensitivity  analysis  is 
performed  to  quantify  the  influence  of  various  model  parameters  and  empirical 
coefficients  on  the  simulation  results.  Predictions  of  the  model  arc  also  examined  for 
varying  water  levels  and  different  initial  profile  slopes. 

5.1  Evaluation  of  K value  with  Laboratory  Experiments 
The  eleven  experiments  described  in  Chapter  4 are  included  in  the  evaluations  of 
the  cross-shore  model  with  the  fixed  average  K value.  The  input  values  of  several 
parameters  in  the  model  are  summarized  below: 

• The  breaking  index,  y = 0.78. 

• The  decay  coefficient,  K , and  stable  wave  height  coefficient,  T , in  the  D’  wave 
model  are  taken  as  0. 1 5 and  0.4,  respectively, 

• The  dissipation  coefficient.  fid , in  the  roller  model  is  set  to  0. 1 . 

• The  limiting  onshore  and  offshore  slopes  for  avalanching  are  0.18  and  0.22, 
respectively. 
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The  deviation  of  the  profile  depth  root-mean-square  errors  (defined  in  Eq.4. 1 ) 
between  the  results  of  the  average  K value  and  the  individual  best-fit  K value  is  chosen  as 
the  model  evaluation  parameter. 


where  Dev  denotes  a non-dimensional  deviation,  RMS  is  the  root-mean-square  error  of 
the  profile  depth,  the  subscript  f and  b represent  the  results  of  the  model  prediction  with 
the  average  and  best-fit  K values,  respectively.  If  the  predictions  with  the  average  K 
value  are  exactly  identical  to  those  of  the  best-fit  K values,  the  deviation  in  Eq.5. 1 is  zero. 
Generally,  the  sign  of  Dev  is  positive  since  the  model  with  best-fit  K values  should  yield 
less  error  than  that  with  fixed  values. 

The  model  results  of  the  root-mean-square  errors  are  presented  in  Table  5-1  for  the 
eleven  experiments.  Figures  5-1  to  5-1 1 show  the  comparisons  of  the  predicted  profiles 
for  these  experiments.  It  is  seen  that  the  predictions  with  the  average  K value  provide  fits 
to  the  measured  profiles  which  are  nearly  as  good  as  those  with  tile  best-fit  K values  for 
each  experiment.  Although  the  best-fit  K values  range  from  -49%  to  72%  of  the  average 
K value,  ten  out  of  eleven  deviations  of  the  toot-mean-square  errors  arc  less  than  10%, 
and  eight  of  the  eleven  are  within  5%. 

The  model  predicts  bar  profiles  in  the  seven  monochromatic  wave  tests,  although 
the  wide  trough  area  in  Case  500  and  501  are  not  produced.  However,  the  model  does 
not  generate  the  double  bars  present  in  the  LWT  experiments.  In  the  random  wave  tests, 
the  predicted  bars  are  not  very  clear  as  in  the  observed  profiles.  The  reason  could  be  that 
the  generated  random  waves  in  the  controlled  laboratory  conditions  are  much  different 
than  those  generated  randomly  in  the  model. 
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5.2  Sensitivity  Analysis  of  Model  Parameters 
Sensitivity  analysis  provides  valuable  information  regarding  the  physical 
implications  of  the  model  parameters  and  their  relative  effects  on  the  results.  In  the 
following,  the  influence  of  the  principal  model  parameters  on  beach  evolution  is 
discussed  by  reference  to  eroded  volume  for  a specific  case  (CRIEP1 6-2). 

5.2.1  Influence  of  Breaking  Index  y 

The  breaking  index  determines  the  offshore  boundary  position  of  the  model. 
Larger  y values  indicate  the  breaking  point  moves  onshore,  and  vice  versa.  Figure  5-12 
illustrates  the  eroded  volume  as  a function  of  time  for  various  values  of  y . A smaller 
value  of  y implies  larger  beach  erosion  due  to  a wider  surf  zone  and  thus  more  material 
movement  from  the  inshore  to  attain  equilibrium. 

5.2.2  Influence  of  I)'  Breaking  Wave  Model  Parameters 

To  quantify  the  importance  of  variations  in  the  wave  height  calculation,  the  wave 
decay  coefficient  K is  varied.  Figure  5-13  shows  the  evolution  of  eroded  volume  for 
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different  values  of  K . A smaller  value  of  K indicates  that  tlie  wave  heights  decay 
slowly  after  breaking  and  thus  propagate  a greater  distance  landward.  From  Eq.  3.31,  a 
smaller  value  of  K gives  a flatter  shape  of  the  equilibrium  beach  profile  and  thus 
requires  redistribution  of  a greater  amount  of  sand  before  equilibrium  is  attained.  On  the 
other  hand,  a smaller  value  of  the  stable  wave  height  coefficient  f results  in  a steeper 
equilibrium  beach  profile  and  a smaller  eroded  volume,  since  less  material  has  to  be 
moved  from  the  inshore  to  attain  equilibrium.  Figure  5-14  illustrates  the  effect  of  varying 
the  stable  wave  height  coefficient  on  the  eroded  volume. 


Figure  5-13.  Effect  of  wave  decay  coefficient  on  eroded  volume. 
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5.2.3  Influence  of  Roller  Model  Coefficient 

The  dissipation  coefficient  pd  in  the  roller  model  determines  the  roller  energy 
dissipation.  A larger  value  implies  steeper  equilibrium  beach  profile  slope  from 
Eq.3.31,  thus  causing  less  volume  erosion  inshore  (Figure  5-15). 


The  equilibrium  energy  dissipation  per  unit  volume  D.  is  a function  of  grain  size 
through  the  profile  parameter  A . A smaller  value  of  D. , occurring  for  finer  grain  size. 
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corresponds  lo  a flatter  equilibrium  beach  profile,  thus  requiring  more  sand  lo  be  moved 
before  equilibrium  is  attained.  In  Ihe  model,  a modification  of  the  A values  (Table  3-2) 
is  made  (Eq.3.34)  to  account  for  the  wave  effects  on  equilibrium  beach  profile.  Figure  5- 
16  illustrates  the  influence  of  varying  equilibrium  energy  dissipation  on  eroded  volume. 
The  three  A values  are:  0.1 19  m10  from  Table  3-2, 0.136  m10  calculated  from  Eq.3.33, 
and  best-fit  value  of  0.1 14  m10  in  the  fitting  process  of  Section  3.2.4. 


Figure  5-16,  Effect  of  profile  scale  factor  on  eroded  volume. 

5.2.5  Influence  of  Wave  Height  and  Period 

The  effect  of  an  increase  in  wave  height  on  beach  erosion  is  obvious,  since  a larger 
wave  hoight  involves  a larger  amount  of  wave  energy  for  the  beach  profile  to  dissipate  in 
a state  of  equilibrium  and  also  displace  the  offshore  boundary  condition  seaward.  The 
eroded  volume  inshore  will  increase  significantly  as  wave  height  increased. 

Figure  5-17  shows  the  evolution  in  time  of  eroded  volume  for  various  wave 
periods.  The  model  is  absolutely  insensitive  to  wave  period.  Although  wave  period 
appears  in  the  roller  model  (Eq.3. 19)  and  the  profile  scale  factor  modification  (Eq.3.33), 
the  effect  of  wave  period  on  profile  evolution  is  extremely  small.  One  reason  may  be  that 
the  linear  wave  theory  for  shallow  water  conditions  is  applied  in  the  numerical  model. 
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However,  il  becomes  very  difficult  lo  program  an  implicil  solution  scheme  if 
intermediate  water  conditions  arc  considered. 


Figure  5-1 7.  Effect  of  wave  period  on  eroded  volume. 

5.2.6  Influence  of  Limiting  Onshore  and  Offshore  Slopes 

These  slopes  arc  introduced  into  the  model  to  determine  slope  avalanching. 


Figures  5-18  and  5-19  illustrate  the  effects  of  varying  limiting  onshore  and  offshore 
slopes,  and  in ^ , on  profile  erosion,  respectively.  It  shows  that  smaller  limiting 


slopes  causes  more  erosion  and  vice  versa.  However,  further  increasing  the  limiting 
offshore  slope,  the  beach  erosion  does  not  change  in  this  case.  The  rational  is  for  a 
smaller  limiting  slope,  more  sediment  need  to  be  redistributed  from  the  dune  and  beach 
area  to  reach  that  slope,  thus  causing  more  erosion. 

53  Effect  of  Varying  Water  Level 

Tests  901  and  91 1 pcrfoimed  in  the  LWT  experiments  were  unique  and  focused  oi 
investigating  the  influence  of  a tide  (variable  water  level)  on  profile  evolution.  A 
sinusoidal  tide  with  amplitude  of  0.45  m and  a period  of  12  hour  was  used  in  Test  91 1 
whereas  the  water  level  was  constant  for  Test  901 . Figure  5-20  shows  the  measured 
profiles  for  the  two  tests  after  40.3  hours  of  wave  action.  Overall,  the  resulting  profile 
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shape  is  quite  similar  in  the  two  cases,  although  the  bar  is  somewhat  larger  for  Test  911. 
This  is  probably  caused  by  the  increased  possibility  for  the  waves  to  erode  higher  up  on 
the  foreshore  due  to  the  water  level  variation. 

This  case  provides  an  opportunity  to  evaluate  model  predictions  of  profile  change 
induced  by  constant  incident  waves  with  a realistic  variation  in  water  level.  The 
predicted  profiles  after  40.3  hours  for  these  two  tests  are  also  shown  in  Figure  5-20.  The 
predicted  profile  with  variable  water  level  has  more  erosion  in  the  foreshore  and  the 
primary  breaking  bar  is  offshore  of  the  bar  position  of  the  constant  water  level  case,  and  a 
small  second  bar  feature  joint  with  the  primary  bar  is  also  produced  by  the  model. 


Figure  5-18.  Effect  of  on-shore  limiting  slope  on  eroded  volume. 


Figure  5-19.  Effect  of  offshore  limiting  slope  on  eroded  volume. 


Figure  5-20.  Comparison  of  observed  and  predicted  profile  evolution  for  variable  water 
5.4  Initial  Beach  Slope  Effect 

Beach  nourishment  projects  arc  usually  constructed  with  an  initial  plane  beach 
slope.  According  to  Dean's  equilibrium  beach  profile  theory,  the  profile  with  an  initial 
slope  steeper  than  the  equilibrium  slope  will  cause  beach  erosion  and  shoreline  recession. 
On  the  other  hand,  an  initial  mild  slope  will  result  in  beach  accretion  and  shoreline 
advancement.  The  model  is  tested  here  to  investigate  the  final  equilibrium  profile  shape 
with  regard  to  different  initial  plane  slopes.  The  wave  conditions  and  sediment  size  are 
held  constant  (H  = 1 . 1 2 m,  T-7.5  s,  h-4.5  m and  d.Q  = 0.27  mm).  The  four  slopes  tested 
are  1/4, 1/10, 1/50,  and  1/100,  respectively.  Figure  5-21  compares  the  predicted  final 
profiles  after  1000  houis  of  constant  wave  action.  Accretionary  profiles  arc  produced 
from  the  mild  initial  slopes  (1/100  and  1/50),  while  erosive  profiles  are  generated  from 
the  steep  initial  slopes  (1/10  and  1/4).  The  final  profiles  have  the  same  slopes  to  the  bar 
crests,  which  indicates  that  the  initial  slope  has  no  effect  on  the  equilibrium  profile  shape, 
but  determines  the  transport  direction.  The  slopes  in  beach  nourishment  projects  usually 
range  from  1:20  to  1:10. 
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Figure  3-21.  Comparisons  of  equilibrium  beach  profile  for  different  initial  beach  slopes. 

5.5  Summary  of  Model  Evaluation 

The  model  results  (n  - 2)  with  the  fixed  average  K value  were  compared  with  those 
with  individual  best-fit  K values  for  the  eleven  large-scale  wave  flume  experiments  used 
for  calibration.  The  deviations  between  the  fixed  K value  and  best-lit  K values  arc  found 
to  be  very  small  (below  10%  in  10  out  of  eleven  cases). 

The  model  is  somewhat  sensitive  to  most  of  the  model  parameters.  One  exception 
is  the  input  wave  period,  to  which  the  model  is  not  sensitive  at  all.  The  varying  water 
level  (such  as  tide,  storm  surge)  on  profile  evolution  can  be  represented  by  the  model.  It 
is  also  found  that  the  initial  beach  slope  has  no  influence  on  the  equilibrium  beach  profile 
shape,  although  it  affects  the  sediment  transport  direction.  This  is  expected  since  this  is  a 
closed  loop  model  and  should  converge  to  a target  profile. 
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Overall,  the  cross-shore  model  with  a nonlinear  transport  relationship  (n  = 2)  has 
the  capacity  to  generate  bar  features  as  well  as  accretionary  profiles.  It  performs 
reasonably  well  for  beach  erosion  in  both  regular  and  irregular  waves.  In  the  next 
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CHAPTER  6 
MODEL  APPLICATION 

In  this  Chapter,  more  large-scale  laboratory  experiment  data  and  field  data  are 
employed  to  examine  the  performance  of  the  cross-shore  model  which  has  been 
developed.  Comparisons  with  predictions  of  SBEACH  and  CROSS  models  are  also 
included. 

6.1  Laboratory  Experiments  Studies 

Measures  of  model  performance  used  for  analysis  of  laboratoiy  experimental  data 
are  based  on  the  statistical  fit  of  the  calculated  profiles  to  the  measured  profile  for  a given 
simulation.  The  residual  parameter.  Res,  used  by  Zheng  and  Dean  (1995),  provides  a 
relative  measure  of  the  difference  between  the  measured  and  calculated  profiles  as 
compared  to  the  amount  of  actual  change  that  has  occurred  over  a given  range  of  the 
profile.  The  residual  parameter  is  given  by: 


where  h is  the  profile  elevation,  Np  is  the  number  of  points  across  the  profile,  the 
subscripts  "p",  “m"  and  “i”  denote  predicted,  measured  and  initial,  respectively.  Smaller 
values  of  Res  correspond  to  a better  fit.  The  minimum  possible  value  of  Res  is  zero, 
which  indicates  a perfect  simulation.  If  no  changes  occurred,  the  value  of  Res  is  1 . Bar 
generation  and  migration  are  of  primary  interest  in  the  laboratory  experimental  results. 

So  just  SBEACH  model  is  used  here  to  compare  with  the  current  model. 
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6.1.1  SUPERTANK  Projecl 

6.1. 1.1  Background 

The  SUPERTANK  project  was  carried  out  in  August  and  September  199!  at  the 
O.H.  Hinsdale  Wave  Research  Laboratory,  Oregon  State  University  (Kraus  and  Smith, 
1994).  SUPERTANK  was  a multi-institutional  effort  to  obtain  detailed  information  on 
beach  profile  response  to  waves  with  special  emphasis  on  surf  and  swash  zone  processes. 
The  project  generated  an  extensive  data  set  consisting  of  wave  transformation,  velocities, 
profile  surveys,  and  sediment  transport  across  the  profile  for  a large  number  of  test  cases. 

The  channel  used  in  the  SUPERTANK  project  was  104  m long,  3.7  m wide,  and 
4.6  m deep.  A digitally  controlled  hydraulic  wave  generator  capable  of  generating 
random  or  monochromatic  waves  was  located  at  one  end  of  the  tank.  The  wave  generator 
faced  a 76-m  beach  profile  constructed  of  600  m3  of  sand.  Sand  for  the  profile  was 
obtained  from  the  Oregon  coast  and  had  a median  grain  size  of  0.22  mm.  Data  collected 
during  the  SUPERTANK  study  were  gathered  by  a dense  array  of  wave  gauges,  current 
meters,  sediment  concentration  profiles,  pressure  gauges,  video  cameras,  and  other 
instruments.  A summary  of  the  data  sets  collected  at  SUPERTANK  is  provided  by  Smith 
and  Kraus  (1995).  During  the  project,  waves  were  run  for  fairly  short  periods  (c.g.  10, 

20, 40,  or  70  min)  in  order  to  avoid  sciching  in  the  tank  and  to  keep  the  amount  of  data 
collected  at  a manageable  level.  Random  waves  with  broad-  and  narrow-banded  spectra 
were  run  with  zero-moment  wave  heights  in  the  range  0.2- 1.0  m and  with  peak  spectral 
periods  in  the  range  3-10  s. 

The  test  conditions  investigated  in  the  present  study  are  summarized  in  Table  6-1 . 
The  developed  model  with  the  recommended  K values  and  other  default  parameters  is 
applied  to  these  tests,  and  the  results  arc  compared  to  those  predicted  by  SBEACH 
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(version  2.0).  The  default  value  for  the  transport  coefficient  in  SBEACH  is  1.75*10* 
m4/N  and  other  input  parameters  arc  also  taken  as  the  suggested  default  values. 

Table  6-1 . Summary  of  the  tests  conditions  in  the  SUPERTANK  project. 


6.1. 1.2  Equilibrium  erosion  (random  waves) 

The  equilibrium  erosion  cases  were  the  first  tests  performed  during  SUPERTANK. 
The  initial  profile  was  a planar  beach  and  was  subjected  to  varying  wave  conditions 
including  both  random  and  monochromatic  waves.  The  first  few  equilibrium  tests 
produced  significant  erosion  at  the  foreshore  and  developed  an  offshore  bar.  During  later 
tests,  the  amount  of  erosion  decreased  as  the  profile  approached  an  equilibrium  shape. 

Figure  6-1  shows  the  initial,  measured  and  calculated  profiles  for  the  first  case, 

PI  A.  The  test  resulted  in  removal  of  material  from  the  upper  profile  and  deposition  in  an 
ofTshore  bar.  Although  the  calculated  results  from  both  models  show  similar  trends 
compared  to  the  measurements,  both  models  overpredict  the  amount  of  erosion  at  the 
foreshore.  Offshore  movement  of  the  bar  is  overestimated.  The  current  model  predicts 
the  bar-trough  system,  whereas  SBEACH  does  not. 

6.1. 1.3  Equilibrium  erosion  (monochromatic  waves) 

Figure  6-2  shows  model  results  for  an  equilibrium  erosion  case  involving 
monochromatic  waves.  For  this  case,  both  models  overestimate  the  amount  of  foreshore 
erosion.  The  location  of  the  bar  is  fairly  well  produced,  although  details  of  the  bar  and 
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1 SBEACH  predicts  better  bar  generation  than  the 


region  of  the  profile.  Wave  breaking  at  the  mound  caused  a significant  amount  of  energy 


mound  erosion.  SBEACH  predicted  the  entire  mound  crest  being  uniformly  lowered. 

Figure  6-4  shows  the  results  for  the  monochromatic  wave  test,  PJC.  The  uniform 
breakpoint  of  the  monochromatic  waves  produced  an  increase  in  elevation  of  the  mound 
crest  corresponding  to  the  formation  of  an  offshore  bar.  A nearshore  bar  also  formed 
from  deposition  of  eroded  foreshore  sediment.  In  contrast  to  the  observed  profile 
changes,  SBEACH  predicts  little  change  in  shape  of  the  offshore  mound  and  no  inshore 
bar.  The  predicted  profile  by  the  current  model  has  large  change  in  the  mound,  with  most 

situated  farther  offshore  than  the  measured  bar. 
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6.1. 1.5  Broad-crcstcd  mound 

This  test  was  similar  to  the  previous  test,  except  that  more  sand  was  added  to  create 
a wider  mound  in  the  ofTshore  to  examine  the  effect  of  mound  width  on  wave  dissipation. 

As  shown  in  Figure  6-5,  the  random  wave  test  case  PKA  produced  a trough  at  the 
offshore  end  of  the  mound,  with  material  being  deposited  on  the  offshore  slope  of  the 
mound.  Little  profile  change  was  observed  elsewhere.  The  SBEACH  simulation  erodes 
the  entire  mound  with  most  of  the  eroded  material  being  deposited  across  the  offshore 
slope  of  the  mound.  As  a result  of  the  lower  mound  height,  the  SBEACH  simulation 
allowed  more  wave  energy  to  reach  the  foreshore,  causing  erosion  to  be  overestimated  at 
the  foreshore.  The  current  model  also  predicts  entire  mound  erosion,  but  approximately 
half  the  material  from  the  mound  is  transported  onshore  to  build  up  the  foreshore  portion 
of  the  profile.  Little  foreshore  erosion  is  predicted  by  the  current  model. 

The  monochromatic  wave  test  (case  PKC)  produced  a trough  on  the  offshore  end  of 
the  broad-crested  mound  and  reduced  the  height  of  the  entire  mound  crest,  as  observed  in 
Figure  6-6.  The  sediment  removed  from  the  mound  was  deposited  in  a bar  at  the  offshore 
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end  of  the  mound  crest  SBE  ACH  produces  slightly  more  lowering  of  the  mound  than 
was  observed,  and  formation  of  a bar  at  the  offshore  end  of  the  mound  is  predicted,  but 
the  observed  trough  is  not  produced.  Similar  to  the  random  wave  case,  the  current  model 
overpredicts  the  mound  erosion  and  onshore  transport  of  the  eroded  material  from  the 
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simulation,  SUPERTANK  I 
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Figure  6-7  presents  the  calculated  residuals  for  the  SUPERTANK  selected  cases. 
From  the  results,  the  predicted  profile  changes  under  both  monochromatic  and  random 
waves  have  some  differences  from  the  laboratory  measurements,  especially  during  the 
offshore  mound  tests.  Some  wave  generated  sediment  transport  mechanisms  are  not 
understood  nor  included  in  both  numerical  models  such  as  wave  plunging,  sediment 
advection  and  diffusion  near  the  breaking  point.  The  residuals  of  current  model  are 
almost  the  same  as  SBEACH  in  P I A,  PGA  and  PJA  cases,  but  have  larger  values  in  PJC, 
PKA  and  PKC  cases  due  to  more  ovcrprcdiction  of  mound  erosion. 
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Figure  6-7.  Residuals  for  SBEACH  and  current  model,  SUPERTANK  experiments. 

6.1.2  LIP  1 1 D Delta  Flume  Experiments 
6.1.2. 1 Background 

In  spring  of  1993  from  April  to  June,  within  the  framework  of  the  European  Large 
Installation  Plan  (LIP),  detailed  measurements  of  hydrodynamics  and  sediment  transport 
in  the  surf  zone  were  undertaken  in  the  Delft  Hydraulics  Delta  Flume.  The  tests  (code 
1 1 D and  thus  denoted  LIP  1 1 D)  were  designed  and  overseen  by  researchers  from  the 
universities  of  Catalunya,  Liverpool,  Padova  and  Thrace.  The  objective  of  the 


experiments  was  to  generate  high-quality  and  high-resolution  data  on  the  hydrodynamics 
and  sediment  transport  dynamics  of  a natural  2DV  beach  under  equilibrium,  erosive  and 
accretive  conditions.  Special  attention  was  paid  to  long  wave  effects  and  near-bottom 
resolution.  A first  overview  of  the  experiments  is  presented  in  Arcilla  ct  al.  (1994).  The 
data  have  been  reported  by  Roelvink  and  Reniers  (1995). 

Two  initial  beach  geometries,  starting  with  a near-equilibrium  Dean-type  beach 
without  (Test  1 ) and  with  the  presence  of  a low  dune  (Test  2),  using  a median  grain 
diameter  of  0.22  mm,  were  used  in  the  project  (Figure  6-8).  The  dune  fool  was  located 
just  above  the  still  water  level  to  investigate  the  effect  of  upper  boundary  condition  on  the 
profile  development  For  each  of  these  two  geometries  the  following  three  wave 
conditions  were  tested:  slightly  erosive,  highly  erosive  and  strongly  accretive.  Narrow- 
banded  random  waves  with  a JONSWAP  spectrum  were  chosen  such  that  the 
combination  of  wave  steepness  at  peak  frequency,  bottom  slope  and  water  level  would 
result  in  a stable,  erosive  and  accretive  beach,  consecutively.  During  each  of  these 
stages,  the  water  level  and  wave  conditions  were  kept  constant.  This  sequence  of 
conditions  could  be  interpreted  as  a schematization  of  the  natural  development  of  a storm, 
represented  by  the  different  sea  states.  The  duration  (number  of  wave  hours)  of  each 
state  was  determined  by  the  lest  duration  and  number  of  possible  wave  hours  per  day,  but 
was  chosen  to  be  at  least  long  enough  to  obtain  accurate  sediment  transport  estimates 
from  the  profile  measurements.  A summary  of  all  test  conditions  is  presented  in  Table  6- 
2,  in  which  incident  wave  heights,  periods,  water  levels  and  durations  arc  described. 
6.1.2.2  Test  1 (without  dune) 

Under  the  first  wave  condition  ( la),  the  initial  profile  changes  into  a mildly  baned 
profile  (Figure  6-9).  The  sediment  transport  within  the  surf  zone  is  predominantly 
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offshore  directed,  while  the  sediment  transport  is  directed  onshore  outside  the  surf  zone. 
During  the  highly  erosive  wave  condition  (lb),  the  bar  region  becomes  more  pronounced 
and  at  the  same  time  the  bar  is  moving  offshore,  as  shown  in  Figure  6-10.  During  the 
accretive  wave  condition  (1c),  the  bar  crest  becomes  more  pronounced,  now  propagating 
in  the  onshore  direction  (Figure  6-11).  The  profile  changes  in  the  trough  and  swash  zone 
are  small  compared  to  the  two  previous  conditions. 

SBEACH  predicts  foreshore  erosion  in  all  three  cases  but  no  bar  is  formed. 
However,  the  sediment  transport  direction  predicted  by  the  current  model  in  the  foreshore 
region  is  predominantly  onshore  for  all  wave  conditions.  The  reason  is  that  the  initial 


profile  slope  is  milder  than  the  anticipated  equilibrium  beach  profile  slope  (Eq.3.31),  thus 
causing  onshore  sediment  movement  according  to  the  EBP  concept. 


Figure  6-8.  Initial  beach  geometry  of  LIP  11D  experiments  (Roelvink  and  Reniers,  1995). 
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6.1.2.3  Test  2 (with  dune) 


The  observed  and  predicted  profile  evolutions  in  the  case  with  a low  dune  are 
presented  in  Figures  6-12  to  6-15.  Starting  with  the  mildly  erosive  wave  conditions,  the 
measured  changes  in  the  bottom  profile  are  similar  to  the  case  without  the  dune.  During 
the  highly  erosive  wave  condition  of  Test  2b,  the  bar  feature  grows  while  moving 
offshore,  increasing  the  water  level  (2c)  results  in  a strong  attack  on  the  low  dune, 
resulting  in  significant  changes  in  the  bottom  profile.  A second  bar  develops  closer  to  the 
shore  whereas  the  outer  bar  decreases  in  amplitude  and  progresses  in  the  onshore 
direction.  Lowering  the  water  level  to  its  original  position  and  using  an  accretive  wave 
condition  (2c)  rebuilds  the  outer  bar  and  diffuses  the  inner  bar. 

As  in  Test  1 , beach  erosion  in  the  foreshore  is  predicted  by  SBEACH  model  in  all 
four  cases  and  a flat  profile  shape  is  generated.  The  current  model  also  predicts  almost 
flat  profiles  in  2a  and  2b  tests,  with  some  dune  erosion  occurring  due  to  the  slope 
adjustment  in  the  model.  There  are  two  bars  reproduced  by  the  current  model  in  case  2e. 
But  the  bar  positions  are  offshore  of  those  observed.  The  model  predictions  in  Test  2c 
arc  similar  to  those  for  Test  lc. 
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Figure  6-16  shows  the  residuals  of  Ihc  model  prediction  for  all  tests  in  LIP  1 ID 
experiments.  The  residuals  of  SBEACH  and  CROSS  are  almost  the  same  except  in  Test 
2a  for  which  the  SBEACH  residual  is  the  smaller  of  the  two.  However,  the  transport  is  in 
opposite  direction  during  the  tests  without  dune. 


6.2  Field  Studies 


Profile  response  at  the  foreshore  and  dune  is  important  in  engineering  applications 
of  numerical  models,  where  model  estimates  of  erosion  are  used  to  determine  the  storm 
effects  on  structures  and  properties  adjacent  to  the  shoreline.  Thus,  measures  of  model 
performance  for  field  cases  focus  on  its  capabilities  to  predict  the  response  of  the  dune 
and  foreshore  after  storms.  These  measures  investigated  in  this  study  include  volume 
change  at  the  foreshore  and  recession  of  a special  contour. 

Volume  change  at  the  foreshore  is  computed  as  the  net  gain  or  loss  in  volume  per 
unit  beach  width  above  some  specified  elevation.  In  this  study,  the  National  Geodetic 
Vertical  Datum  (NGVD)  is  used  as  the  reference  elevation.  Contour  recession  provides 
uuumiauuii  on  now  well  a model  predicts  the  erosion  of  a dune  or  berm  and  the  loss  of 
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reference  elevation  (NGVD)  and  the  dune  crest  is  selected  for  each  field  case.  Recession 

Field  cases  selected  for  this  study  arc  based  on  the  availability  of  data  records  of 
beach  profile  surveys,  wave  and  water  level  information,  and  sediment  grain  size 
information.  The  most  recognized  and  applied  field  data  set  in  the  literature  is  Ocean 
City,  Maryland  (a  scries  of  storms  from  October  1991-January  1992)  (Stauble  el  al„ 
1993).  The  profile  changes  after  Hurricane  Eloise  1975  and  Hurricane  Opal  1995  in  the 
Panhandle  Coast  of  Florida  are  also  investigated. 

6.2.1  Ocean  City,  Maryland 

6.2.1. 1 Background 

Ocean  City,  MD,  is  located  on  Fenwick  Island,  a north-south  oriented  sandy  barrier 
Island  of  the  central  Delaware-Maryland- Virginia  coast.  The  island  extends  from  Indian 
River  Inlet,  Delaware,  on  the  north  to  Ocean  City  Inlet  on  the  south  and  is  backed  by  Isle 
of  Wight  and  Assawoman  Bays.  The  Ocean  City  beach  was  nourished  by  the  state  of 
Maryland  in  1988,  and  the  Federal  Government  in  1990  and  1991  to  protect  the  city 
against  storm  damage  (Stauble  ct  al„  1993).  The  project  layout  is  shown  in  Figure  6-17. 

and  3.  The  representative  grain  sizes  from  Borrow  Areas  2 and  3 arc  0.25  mm  and  0.35 
mm,  respectively.  The  fill  constructed  by  the  Stale  of  Maryland  consisted  of 
approximately  2.7  million  cubic  yards  (2.1  million  cubic  meters)  of  material  placed  along 
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million  cubic  yards  (2.9  million  cubic  meters). 


: 6-17.  Tile  Ocean  City  be 


I (St  i 


al„  1993). 
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From  late  October  1991  to  early  January  1992  after  the  project,  the  beach  at  Ocean 
City  was  impacted  by  a series  of  severe  storms.  The  series  of  storms,  which  included  the 
October  30, 1991  "Halloween  storm",  November  11, 1991  storm,  and  January  4, 1992 
storm,  caused  extensive  erosion  of  the  constructed  foreshore  and  dune  at  Ocean  City.  As 
part  of  the  project  monitoring  effort,  beach  profile  data  at  ten  locations  along  the 
shoreline  were  collected  which  enabled  the  analysis  of  the  erosional  effects  of  the  series 
of  storms.  The  wave  height,  wave  period  and  storm  surge  time  history  were  measured  by 
two  gages  located  directly  offshore  of  Ocean  City  in  a water  depth  of  10  m.  Water  depths 
and  profile  elevations  used  are  referenced  to  NGVD,  which  lies  0.02  m below  mean 
water  level  for  Ocean  City- 

Six  survey  lines  located  at  37th,  45th,  56th,  63rd,  74th,  and  103rd  Streets  have  both 
pre-storm  and  post-storm  measured  profiles  available  and  are  simulated  in  this  study  for 
all  three  major  storms.  In  the  analysis  of  the  profile  changes,  the  net  volume  gained  or 
lost  between  surveys  are  quite  different  from  zero,  which  is  considered  to  be  due  to 
gradients  in  longshore  sediment  transport.  To  remove  this  effect,  each  post-storm  profile 
is  adjusted  by  shifting  the  profile  horizontally  a distance  Ay  to  yield  zero  net  volume 
change  (Zheng,  1996;  Zheng  and  Dean,  1997).  The  value  of  Ay  is  calculated  by 

4V  = — /(*_  (6.2) 

where  subscripts  mb  and  ma  denotes  profile  elevation  measured  before  and  after  storms, 
respectively,  y0  and  y,  are  the  offshore  distance  coordinates  at  the  baseline  and  offshore 
profile  change  limit,  respectively,  and  is  the  total  elevation  of  the  active  post-storm 
profile.  The  sign  of  Ay  is  positive  for  a seaward  translation. 
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Based  on  results  of  sediment  analyses  performed  during  the  monitoring  program,  a 
median  sediment  grain  size  of  0.35  mm  is  selected  as  a representative  value  for  all  beach 
profiles  in  the  numerical  simulation, 

6.2.1, 2 Halloween  storm 

The  northeast  storm  of  October  29  - November  2, 1991  (the  “Halloween  Storm") 
produced  high  waves  and  water  levels  at  Ocean  City,  with  a peak  significant  wave  height 
of  about  3 m and  a peak  water  elevation  of  about  1 .5  m NGVD.  Figure  6- 1 8 shows  the 
measured  significant  wave  height,  peak  spectral  wave  period  and  storm  surge  time 
histories  for  this  storm.  Pre-storm  profiles  were  collected  in  June  1991  and  post-storm 
surveys  were  performed  shortly  after  the  storm.  The  June  surveys  were  assumed  to 
adequately  represent  the  conditions  of  the  upper  beach  profiles  prior  to  the  storm,  as  no 
other  significant  erosion  events  occurred  between  June  and  October.  Table  6-3  presents 
the  net  volume  change  per  unit  beach  width  calculated  between  surveys  for  each  profile 
and  the  corresponding  beach  shifting  adjustment  calculated  from  Eq.6.2.  Five  profiles 
exhibited  a net  loss  of  material  while  one  profile  (103rd  Street)  gained  material. 

The  model  predictions  with  default  input  parameters  for  SBEACH.  CROSS  and 
current  model  ate  presented  in  Appendix  B for  all  six  profiles.  Figures  6- 19  and  6-20 
show  one  example  of  the  comparison  of  the  calculated  results  for  the  profile  at  63rd  street 
with  the  original  and  shifted  profile,  respectively.  Good  agreement  is  obtained  above  the 
NGVD  datum.  Below  datum,  the  measured  profile  exhibits  a convex  shape  with  a 
nearshore  step,  whereas  the  calculated  profiles  by  all  three  models  show  more  concave 
shape  with  a bar  situated  offshore  by  SBEACH  and  the  current  model.  The  measured 
beach  step  may  be  attributed  to  beach  recovery  prior  to  the  surveys  during  the  waning 
storm  conditions,  an  effort  not  included  by  the  models. 
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Figure  6-19.  Comparison  of  calculated  and  original  measured  profile  of  Ihe  Halloween 
slorm  for  the  profile  at  63rd  street.  Ocean  City,  MD. 
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Figures  6-21  to  6-23  compare  volume  changes  above  NGVD,  recessions  of +1 .5  m 
contour,  and  residuals  for  the  measured  and  simulated  profiles  for  each  profile  (street), 
respectively.  In  the  figures,  "without  adjustment"  means  data  given  by  the  original 
measured  post-storm  profiles,  while  "with  adjustment"  means  data  resulting  from 
horizontally  shifted  post-storm  profiles.  The  calculated  volume  changes  and  contour 
recession  by  all  models  arc  closer  to  the  measured  results  with  adjustment.  On  average, 
the  current  model  results  are  similar  to  those  by  SBEACH.  Some  differences  between 
SBEACH  and  the  current  model  may  be  due  to  the  current  model  not  including  sediment 
transport  in  swash  zone  and  overwash,  which  are  incorporated  in  SBEACH  model.  The 
residuals  by  CROSS  are  smaller  than  the  other  two  models  which  may  be  interpreted  that 
CROSS  has  the  belter  overall  profile  prediction.  However,  after  examining  the  calculated 
profiles,  it  is  believed  that  the  major  differences  are  caused  by  the  offshore  profile  shape 
predictions.  CROSS  cannot  produce  bars  and  generates  a fiat  offshore  profile,  which 
occurs  in  the  measured  profiles  in  coincident  due  to  beach  recovery.  The  largest  residual 
discrepancy  is  at  the  56th  street,  where  CROSS  predicts  minor  profile  changes  while 
SBEACH  and  the  current  model  have  beach  erosion  at  the  pre-storm  nearshore  step  area 
and  deposition  offshore. 

6.2.1 3 Novcmber-Januan  storms 

The  November-January  (NJ)  storms  include  the  November  11, 1991  storm  and  the 
January  4, 1992  storm,  with  the  measured  significant  wave  height,  peak  wave  period  and 
storm  surge  level  shown  in  Figures  6-24  and  6-25,  respectively.  In  this  case,  these  two 
storms  arc  combined  together  in  the  simulation  because  no  beach  profile  measurement 
data  are  available  between  these  two  storms. 
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Ocean 


Figure  6-24.  Significant  wave  heighl,  peak  wave  period  and  waier  level  lime  history  for 
the  November  1991  storm.  Ocean  City,  MD  (Stauble  et  al„  1993). 
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As  observed  in  the  figures,  high  waves  and  wafer  levels  lasted  approximately  3 
days  for  both  storms.  The  November  storm  had  a peak  significant  wave  height  of  about  3 
m and  water  elevation  of  about  1 .2  m NGVD.  The  January  storm  is  much  more  severe, 
with  a peak  wave  height  and  water  level  of  4 m and  2 m NGVD,  respectively.  The 
meteorological  characteristics  of  the  January  1992  storm  contrast  significantly  with  those 
of  the  Halloween  Storm.  Whereas  the  Halloween  Storm  was  a very  large,  slow-moving 
maritime  winter  cyclone,  the  January  storm  was  a small,  rapidly  developing,  fast-moving 
event  (Staublcctal.,  1993). 

The  beach  profiles  measured  after  the  Halloween  storm  served  as  pre-storm  profiles 
for  the  NJ  storm  set.  Beach  surveys  performed  shortly  after  the  January  storm  provided 
post-storm  profiles  for  comparison  with  model  results.  Table  6-4  presents  the  net  volume 
change  per  unit  beach  width  calculated  between  surveys  for  each  profile  and  the 
corresponding  beach  shifting  adjustment  after  the  NJ  stoim.  A large  gain  in  material  was 
observed  on  profiles  at  37th  and  103rd  street,  whereas  volume  was  conserved  reasonably 
well  on  the  other  profiles. 

Measured  beach  response  to  the  NJ  storm  scries  shows  that  four  of  the  six  profiles 
experienced  overwash  of  the  dune,  and  SBEACH  produces  dune  overwash  in  each  of  the 
four  cases,  whereas  CROSS  and  the  current  model  do  not.  Figures  6-26  and  6-27 
illustrate  the  comparison  of  simulated  results  for  the  profile  at  63rd  street  with  original 
and  shifted  profiles,  respectively.  It  is  noted  that  although  substantial  erosion  occurred  at 
the  dune,  erosion  on  the  lower  foreshore  profile  was  minimal.  A berm  above  NGVD  was 
demonstrated  similar 


e.  Profiles  at  45th,  74th,  and  103rd  street  i 


Tabic  6-4.  Measured  volume  changes  during  the  NJ  sic 
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Volume  changes  above  NGVD,  recessions  of +1.5  m contour,  and  residuals  of 
profile  changes  are  presented  respectively  in  Figures  6-28, 6-29  and  6-30  for  the  NJ 

SBEACH  results  and  the  differences  could  be  improved  by  incorporating  swash  zone 

subaerial  at  the  56th  street  profile  in  both  observed  and  predicted  results  by  SBEACH  and 
CROSS,  which  is  indicated  by  volume  change  and  contour  recession.  The  current  model 
predicts  some  onshore  sediment  movement  from  the  pre-storm  profile  near  NGVD  and 

underpredicted  by  all  three  models  at  the  63rd  and  74th  streets,  which  are  located  in  the 
area  of  the  crosional  hotspot  (Wise  et  al.  1996).  The  better  residual  results  from  CROSS 
are  believed  to  be  the  same  reason  as  mentioned  for  the  Halloween  storm. 

6.2.2  Panhandle  Coast  of  Florida 

Several  severe  storms  have  struck  the  Panhandle  Coast  ofFIorida  (Figure  6-31)  the 
past  century.  Here  beach  erosion  caused  by  Hurricanes  Eloise  and  Opal  are  examined. 
6.2.2.1  Hurricane  Eloise 

the  Dominican  Republic.  Cuba,  and  adjacent  waters,  and  regaining  its  hurricane  force  on 
the  morning  of  September  22, 1975  in  the  central  Gulf  of  Mexico  about  300  nautical 
miles  south  of  New  Orleans,  struck  the  Florida  northwest  coast  in  the  early  morning  of 
September  23, 1975.  Its  center  crossed  the  shoreline  about  midway  between  Fort  Walton 
Beach  in  Okaloosa  County  and  Panama  City  Beach  in  Bay  County  at  a forward  speed  of 
better  Ulan  20  knots  towards  the  northeast.  Hurricane  Eloise  was  classified  3 on  a scale 
of  5 by  the  National  Hurricane  Center. 
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Hydraulic  conditions.  Hurricane  Eloise  was  a rapidly  moving  storm,  lasting  less 
than  20  hours.  There  are  no  open  coast  storm  surge  measurements  available.  Numerical 
estimates  ranged  from  8.4  feet  (2.56  m)  at  the  western  end  of  Walton  County  to  9.6  feet 
(2.93  m)  at  the  eastern  end  of  the  county  (Kriebel,  1986).  Figure  6-32  shows  the 
predicted  storm  surge  hydrograph  near  the  western  end  of  Walton  County  used  by 
Kriebel  (1986)  to  verify  his  dune  erosion  model.  The  storm  surge  at  other  sites  was 
determined  by  multiplying  a constant  over  the  surge  hydrograph  in  Figure  6-32  for  each 
profile.  A representative  significant  breaking  wave  height  of  1 2 feet  (3.66  m)  and  peak 
period  of  1 1 seconds  were  used  in  Kriebel  (1986). 


Figure  6-31.  Panhandle  Coast  ofFlorida  (Leadon  ct  al.,  1998). 
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Profile  response.  Under  the  Florida  Coastal  Construction  Setback  Line  Program, 
the  coastal  areas  primarily  affected  by  Hurricane  Eloise  in  Bay,  Walton  and  Okaloosa 
Counties  were  surveyed  in  October  1973,  with  beach-dune  profiles  at  approximately 
every  1000  feet,  and  offshore  soundings  every  3000  feet.  Post-storm  surveys  were 
conducted  within  2-4  weeks  after  the  storm  in  October  1975  in  Walton  County. 

The  entire  Walton  County  had  suffered  severe  bcach-dunc  erosion  during  the 
storm.  However,  the  measured  profiles  showed  Oral  there  was  great  variability  in  pre- 
storm  profile  forms  and  in  storm  erosion  response  of  individual  profiles.  In  areas  where 
the  beach  was  wide  with  gentle  slope  and  gradually  merged  into  the  dunes,  there  was 
little  erosion;  whereas  in  the  neighboring  areas  with  steeper  beach  slope  and  rapid  dune 
rise,  severe  erosion  took  place.  Such  instances  indicate  that  gradual  wave  energy 
dissipation  reduces  the  erosion  potential  on  the  beach.  Also  over  20  profiles  showed  a 
net  gain  of  sand  between  the  two  available  surveys.  Such  variability  in  profile  response 
con  be  attributed  in  part  to  three  sources  of  uncertainty;  (1)  pre-storm  profile  modification 
between  October  1973  and  Hurricane  Eloise,  (2)  localized  longshore  transport  effects 
during  the  storm,  causing  areas  of  local  accretion  and  accentuated  erosion,  and  (3)  non- 
uniform  beach  recovery  between  the  storm  and  the  post-storm  surveys.  A sand  ridge  near 
and  above  the  mean  sea  level  on  the  post-storm  profiles  in  Walton  County  was  found  to 
be  about  1. 5-2.5  feet  (0.46-0.76  m)  high  and  50  feet  (15  m)  wide.  On  average  50  fl3/ft 
(4.7rtr/m)  had  returned  to  the  beach  face  prior  to  the  post-storm  survey  (Chiu,  1977).  No 
offshore  bar  was  formed  on  the  post-storm  profiles  which  may  be  because  the  storm  was 
so  fast-moving  that  there  was  not  enough  time  to  form  the  bar  or  a barely  formed  bar  was 
eliminated  by  waves  between  the  storm  and  the  survey. 
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Numerical  simulation.  A reference  profile,  line  R-41,  which  has  been  used  in 
other  dune  erosion  verification  studies  (Hughes  and  Chiu  1981;  Velliga  1983;  Kriebel 
1986),  is  chosen  to  investigate  the  model  performance  on  Hurricane  Eloise.  Due  to  lack 
of  offshore  profile  measurement,  the  pre-storm  profile  is  constructed  as  in  Kriebel 
(1986):  above  NGVD,  the  actual  measurement  is  applied;  offshore,  a Ay11’  profile  is 
assumed.  The  scaling  parameter  A,  based  on  the  effective  grain  size  of 0.262  mm,  is 
0.184  ft10  (0.123  m1'3).  The  hydraulic  conditions  of  Kriebel  (1986)  are  also  adopted. 

The  pre-storm  and  post-storm  profiles  are  shown  in  Figure  6-33.  All  three  models 
uiidcrcsuinaic  die  dune  erosion  substantially.  The  above-mentioned  first  two  factors  in 


6.2.2.2  Hurricane  Opal 

Background.  Hurricane  Opal,  which  struck  the  Panhandle  Coast  of  Florida  on 
October  4, 1 995,  was  one  of  the  most  severe  hurricanes  to  impact  Florida  during  the  past 
century.  Opal  reached  category  4 status  as  it  approached  landfall  on  October  4, 1995 


wilh  maximum  sustained  surface  winds  reaching  ISO  mph.  Hurricane  Opal  caused  more 
structural  damage  along  the  Florida  coast  than  all  hurricanes  and  tropical  storms 
combined  since  1975.  Severe  damage  occurred  across  the  entire  Panhandle  region  with 
coastal  flooding  extending  along  the  central  Florida  Gulf  coast  region  southward  all  the 
way  through  Key  West. 

Hydraulic  conditions.  There  were  several  NOAA  stations  that  recorded  data  as 
Opal  passed  through  the  Gulf  and  onto  land  (Figure  6-34).  Wave  data  collected  at  the 
offshore  buoy  station  42036,  which  is  nearest  to  the  Florida  coast,  is  shown  in  Figure  6- 
35.  The  data  at  this  station  shows  the  dramatic  increase  in  wave  height  and  dominant 
wave  period  as  Opal  traveled  across  the  Gulf,  even  though  Opal’s  track  was  quite  a 
significant  distance  west  of  this  station. 


Panama  City  Beach  Pier  (near  Monument  line  R-40)  is  shown  in  Figure  6-36.  The  peak 
storm  surge  elevation  is  shown  at  +2.53  m above  the  NGVD  datum.  The  tide  gage 
collected  this  data  near  the  seaward  end  of  the  pier,  several  hundred  feet  seaward  of  the 
shoreline.  A storm  surge  model  estimated  a peak  storm  surge  at  the  shoreline  to  be  from 
+3.7  to  +4.9  m (NGVD)  (Lcadon  et  al.,  1 998),  which  is  consistent  with  a number  of  high 
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water  marks  surveyed  in  the  Panhandle  by  FDEP  staff  and  other  surveyors  following 


Figure  6-36.  Storm  surge  hydrograph  for  Hurricane  Opal  at  Panama  City  Beach  (Leadon 
etal.,  1998). 

waves  severely  eroded  the  beach  and  dune  system  throughout  the  Panhandle  Coast. 
Beach  erosion  resulted  in  a significantly  lower  beach  and  a shoreline  that  retreated 


areas  with  lower  dune  topography  experienced  dune  overtopping  in  which  sand  material 
was  transported  inland  as  overwash. 

Beach  and  offshore  profile  data  were  collected  by  the  Bureau  of  Beaches  and 
Coastal  Systems,  FDEP,  in  March  of  1 995  in  Bay  County  and  in  June  of  1 995  in  Walton 
County.  The  post-storm  beach  and  dune  profiles  data  were  collected  in  October  1995,  in 
many  cases  within  4-5  days  of  Opal's  October  4th  landfall.  In  addition,  the  Mobile 
District  of  the  U.S.  Army  Corps  of  Engineers  obtained  a high-resolution  bathymetric 
survey  of  the  offshore  region  in  the  vicinity  of  their  planned  beach  restoration  project  at 
Panama  City  Beach  (unfortunately  these  data  arc  unavailable).  Typical  examples  of 
profile  comparisons  for  Bay  County  which  include  the  offshore  survey  are  presented  in 
Figure  6-37  (original  Figures  43  and  44  in  Leadon  el  al.  1998).  The  detailed  analysis  of 
beach  and  dune  erosion  and  structural  damage  along  the  Panhandle  Coast  after  Hurricane 
Opal  are  given  by  Leadon  et  al.  (1998). 

Numerical  simulation.  Several  profiles  from  Panama  City  Beach  in  Bay  County 
are  selected  to  compare  the  numerical  model  predictions  with  measured  profiles.  The 
sand  grain  size  is  taken  as  0.26  mm.  The  WIS  data  at  Station  182  (Figure  6-38)  is 
coupled  with  the  storm  surge  information  in  Figure  6-36  as  the  input  hydraulic 
conditions.  Figures  6-39  and  6-40  show  two  examples  of  the  results  from  profile  R-21T 
and  R-41,  respectively.  As  in  the  Hurricane  Eloise  case,  all  three  models  underestimate 
the  dune  erosion  substantially.  The  outer  bars  on  the  pre-storm  profiles  are  not  totally 
eroded  by  SBEACH  and  the  current  model.  However,  in  the  observed  profile  shown  in 
Figure  6-37  for  R-21T,  the  original  bar  was  totally  smoothed  out  by  the  storm. 
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nodelsSBEACH  and  CROSS.  Quantities  chos 


model  performance  include  residual  of  profile  change,  volume  change  above  selected 
elevation,  and  recession  of  a selected  contour. 

In  the  laboratory  experiments  studies,  because  the  focus  is  on  the  bar  profile 
generation  and  migration  under  different  wave  conditions,  and  CROSS  model  could  not 
represent  bar  profiles,  only  comparison  with  SBEACH  and  residuals  of  profile  change  arc 
investigated.  For  the  simulation  of  storm  erosion  at  Ocean  City,  MD,  an  adjustment  is 
made  by  shifting  the  entire  profile  horizontally  a distance  to  yield  zero  net  volume  change 
to  approximate  removal  of  longshore  sediment  transport  effects.  This  adjustment  may 
not  be  appropriate  for  the  upper  beach  profile  and  dune  region  because  the  longshore 
processes  primarily  affect  the  subaqueous  profiles.  From  the  results,  it  appears  that  the 
residuals  are  less  affected  by  the  shifting  adjustment,  while  the  eroded  volumes  and 
contour  recessions  arc  affected  more  significantly  by  the  shift.  For  Hurricanes  Eloise  and 
Opal,  there  is  no  profile  adjustment  due  to  lack  of  offshore  profile  measurements. 

The  results  show  that  the  current  model  with  default  parameters  provides  fairly 
reasonable  agreement  in  most  cases  with  the  measurements  and  model  results  of 
SBEACH.  Both  models  could  produce  bar  features  offshore,  although  the  bar  crest  is 
more  pronounced  in  the  current  model  than  SBEACH.  The  differences  between  CROSS 
and  the  current  model  and  SBEACH  are  believed  to  be  mainly  due  to  different  wave 
energy  dissipation  representations.  CROSS  uses  the  simple,  spilling  breaker  assumption 
inside  the  surf  zone,  while  the  current  model  and  SBEACH  incorporate  more  physically 
sound  dissipation  mechanisms.  The  discrepancy  of  volume  losses  of  all  three  models  on 
field  data  is  believed  to  be  somewhat  mainly  due  to  longshore  transport  effects  during 


CHAPTER  7 

SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 
7.1  Summary 

A numerical  model  is  developed  to  investigate  the  use  of  a roller  model  for  bar 
generation  in  the  nearshore  region  and  to  better  represent  the  lime  history  of  cross-shore 
profile  evolution.  The  cross-shore  sediment  transport  relationship  in  the  model  is  based 
on  equilibrium  beach  profile  concepts.  The  wave  transformation  model  by  Dally  ct  al. 
(1984)  and  the  roller  model  by  Dally  and  Brown  (1995)  are  applied  to  predict  the  wave 
height  and  roller  energy  dissipation  inside  the  surf  zone,  respectively.  Profile  evolution 
under  random  waves  is  simulated  by  the  wave-by-wave  approach  which  does  not  require 
assumptions  of  wave  height  distribution  inside  the  surf  zone.  In  order  to  represent  bar 
profiles,  a grid  system  is  used  that  can  include  representation  of  bar  features  readily.  An 
implicit,  double  sweep  scheme  is  adopted  to  solve  the  sediment  transport  and  continuity 
equations  numerically.  An  effort  is  made  to  account  for  the  wave  effects  on  equilibrium 
beach  profiles  using  the  fall  velocity  parameter  (H/wT).  The  two  unknowns  in  the  cross- 
shore sediment  relationship  arc  the  transport  coefficient,  K,  and  the  exponent,  n,  in 

Q - K{D  - D.  )]£>  - 17.  j"  . Up  to  eleven  large-scale  wave  flume  laboratory  experiments 
from  three  different  facilities  arc  incorporated  to  calibrate  and  compare  the  sediment 
transport  relationships  (K  and  n = 1 , 2 and  3).  The  best-fit  transport  coefficients  arc 
determined  as  the  value  yielding  the  overall  least  square  error  of  profile  changes  in  each 
case.  Based  on  the  calibration  results,  it  is  found  that  the  non-linear  transport  relationship 
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with  n = 2 provides  the  best  representation  of  profile  evolution  for  monochromatic  waves, 
a result  reinforced  in  investigating  eroded  volume  evolution  with  time.  For  the  four 
irregular  wave  tests,  n - I provides  the  best  fit. 

The  performance  of  the  model  is  evaluated  by  the  difference  between  the  errors  of 
profile  change  predicted  by  the  model  with  fixed  and  best-fit  K values  using  the  same 
laboratory  experimental  data  as  in  the  model  calibration.  The  predictions  of  the  model 
with  fixed  and  best-fit  K values  are  quite  close  for  the  eleven  experiments.  Sensitivity 
analyses  are  also  conducted  to  investigate  the  effects  of  several  model  input  parameters 
on  the  model  performance.  The  model  is  then  applied  to  other  large-scale  wave  flume 
experiments  (SUPERTANK  and  UP  1 ID)  and  field  data  (Halloween  storm  1991, 
November  1991  and  January  1992  storm  at  Ocean  City,  Maryland;  Humcane  Eloise  and 
Opal  at  the  Panhandle  Coast  of  Florida).  The  model  predictions  are  compared  with  two 
other  models  (SBEACH  and  CROSS)  with  default  input  parameters  for  all  three  models. 
In  the  Ocean  City  simulation,  in  order  to  reduce  longshore  transport  effects,  an 
adjustment  is  made  by  shifting  the  entire  active  profile  horizontally  a distance  to  yield  a 
zero  net  volume  change  for  each  post-storm  profile.  In  the  Hurricane  Eloise  and  Opal 
cases,  there  is  inadequate  offshore  profile  measurements  to  conduct  such  an  adjustment. 
Three  different  measures  are  proposed  to  evaluate  the  agreement  between  measurements 
and  predictions:  residual  parameter  for  profile  change,  volume  change  and  contour 
recession  above  NGVD. 

7.2  Conclusions 

The  present  numerical  model  is  able  to  simulate  beach  profile  evolution  under  both 
monochromatic  and  random  wave  conditions.  Based  on  the  investigation  conducted  and 
results  developed  herein,  the  following  conclusions  are  offered: 
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(1)  The  inclusion  of  a surface  roller  model  with  a realistic  breaking  wave  model 
provides  a rational  basis  for  bar  generation.  The  roller  model  captures  the  lag  between 
initial  wave  breaking  and  maximum  wave  energy  dissipation.  A new  equilibrium  beach 
profile  shape  is  developed,  which  can  represent  bar  profiles. 

(2)  Modification  of  the  profile  scale  factor  to  include  a dependency  on  the  fall 
velocity  parameter  (H/wT)  accounts  for  the  effects  of  wave  characteristics  on  beach 
profile  shape  and  provides  an  explanation  for  seasonal  shoreline  changes  in  nature. 

(3)  The  exponent,  n = 2,  in  the  sediment  transport  equation  provides  better 
agreement  to  the  laboratory  data  than  values  used  in  previous  closed  loop  models  (n  - I 
and  3)  under  monochromatic  waves.  The  determined  transport  relationship  better 
represents  time  scales  of  profile  evolution  under  different  initial  profile  and  wave 
conditions.  The  transport  coefficient  K of  0.93x  I O'8  m6s/N2  is  recommended  for  model 
application. 

(4)  The  transport  relationship  with  n = 1 provides  the  best-fit  to  random  wave  tests. 
The  reason  probably  is  that  the  generated  waves  using  the  wave-by-wave  approach  are 
different  than  the  actual  wave  series  generated  in  the  wave  flumes,  although  this  issue 
requires  further  investigation. 

(5)  The  model  is  somewhat  sensitive  to  most  of  the  model  input  parameters  with 
the  exception  of  wave  period.  The  model  insensitivity  to  wave  period  may  be  due  to  the 
shallow  water  wave  environment  in  which  this  model  is  applied.  It  is  also  found  that  the 
initial  beach  slope  has  no  influence  on  the  final  equilibrium  beach  profile  shape,  although 
it  affects  the  sediment  transport  direction  and  the  time  required  for  equilibration. 


(6)  The  units  of  the  sediment  transport  coefficient  K,  when  interpreted  by  Froude 
scaling,  has  units  of  (length)1'2.  Thus,  when  applying  this  model  to  field  situations,  it  is 
advisable  to  calibrate  the  K value  if  possible  or  to  consider  scaling  based  on  consideration 
of  field  and  model  scales. 

(7)  Available  field  data  contain  an  unknown  amount  of  recovery  due  to  onshore 
sediment  transport  occurring  before  the  post-storm  surveys  and  some  net  volume  changes 
along  the  profiles  due  to  longshore  sediment  transport  gradients. 

7.3  Recommendations  for  Future  Study 

Several  features  of  the  present  work  merit  further  consideration  and  can  be  readily 
invesugatea  or  incorporated  into  the  model  for  investigation.  The  profile  changes 
between  the  wave  set-up  limit  to  the  runup  limit  (the  swash  zone)  in  the  model  are  treated 
to  adjust  the  slope  once  the  limiting  slope  is  reached.  To  better  represent  the  profile 
changes  in  this  region  and  its  subsequent  effects  on  surf  zone  profile  changes,  the  detailed 
swash  mechanisms  and  the  corresponding  sediment  transport  need  to  be  represented. 

By  incorporating  surface  roller  into  the  wave  energy  dissipation,  a bar-trough 
profile  can  be  generated.  However,  the  detailed  profile  evolution  near  the  bar-trough  area 
is  not  well  simulated.  On  observed  profiles,  shoreward  slope  of  the  bar  is  much  steeper, 
and  the  bar  is  much  wider  than  predicted  by  the  model.  In  some  cases,  multiple  bars  arc 
formed  inside  the  surf  zone.  Some  adjustments  on  sediment  transport  near  the  wave 
breaking  point  may  provide  improvement.  The  convergence  between  offshore  transport 
inside  the  surf  zone  and  on-shore  transport  outside  the  surf  zone  due  to  wave  asymmetry 
may  contribute  to  bar  growth.  Lack  of  transport  representation  outside  of  the  surf  zone 
also  limits  the  model  capability  to  simulate  on-shore  bar  migration  and  beach  recovery  in 


summer  wave  conditions. 


Another  model  improvement  should  be  to  belter  represent  the  wave  period  effects 
on  profile  evolution.  From  the  investigation  on  measured  profiles  under  different  period 
waves  in  Chapter  1 . it  is  clear  that  wave  periods  have  an  influence  on  cross-shore 
sediment  transport,  especially  long  period  waves  (as  in  Figure  1-7).  This  effect  needs 
further  investigation  and  incorporation  into  the  model,  although  the  effect  of  wave 
parameters  (HAvT)  on  the  profile  scale  parameters.  A,  represents  one  contribution  to  this 

The  sediment  transport  and  resulting  profile  change  under  random  waves  arc 
treated  by  the  wavc-by-wave  approach,  which  has  the  advantage  of  no  wave  height 
distribution  assumption  requirements  inside  the  surf  zone.  However,  wave  interaction 
effects  on  sediment  transport  are  not  taken  into  account.  Some  studies  have  shown  that 
long  waves  (infragravity  waves)  generated  by  wave  interactions  can  have  a significant 
influence  on  sediment  transport.  Further  effort  is  required  to  determine  the  best  approach 
of  including  this  effect  in  closed  loop,  energy  dissipation  based  models. 

The  representations  of  cross-shore  sediment  transport  and  profile  evolution  under 
two  other  commonly  met  conditions,  overwash  caused  by  a severe  storm  and  beach- 
structure  interaction  in  front  of  shore  parallel  structures,  are  requirements  for  a model  to 
be  labeled  as  a complete  and  "good"  beach  profile  evolution  model.  Although  these  two 
issues  have  been  included  in  the  SBEACH  model,  they  were  treated  in  ad  hoc  fashions. 

A rational  understanding  of  these  processes  and  incorporating  them  into  the  current 
model  is  another  main  focus  in  further  model  development. 


FORTRAN  PROGRAM  CODE  LISTING 
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REAtM.'ITOirrm  !TIME  TO  OUTPUT  PREDICTED  RESULTS 


READ(V)DIST(I>.DPT(I)  ! INITIAL  PROFILE 

DISTID-DISTIO-LINIT 

OPTOWDPTIIteSWLlnJNIT 


177 


IF(.NOT.ALLOCATED(DISn))  ALL0CATE(DISn(IIMAX*1 )) 
c'fL^PRO^NTEIMlD^^X^G^'^^A^^P^'^n^ppTi 

NGRID-IIMAX  ^ 

!WRITE(6J6KDISn(lyUNIT.-Dm(IVUNrr.l=l.NGRID) 
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DPTCID-DPTCID'O. 


CALLCOEFF(A,B.C.D.AS,BS.CS1DS.DPrC.OO,NGRID*l.WAVH.DlSP.AROLUDSTARP> 

CAU.SSWEEP(E.F.kspS.Dm'DIWBC.lSNGfu™|I!>,l'IBC'IECI 


• EXAMINE  WHETHER  OH<l 

DHMAX-0.0 

IF(ABS(DH(I)).GT.DHMAX)  DHMAX-ABS<DH|I» 
IF(DHMAX.GE.I.O)  STOP  13H*I,  REDUCE  TIME  STEPP 
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if<™ceq.tout<jtd™en 

WRITR6,36XDIST1(1)/UNIT.-DTO(IVU1 
WRITE(7.36KD1STI(I),Q(I).I- 1 .NGRID) 
WRITE|8,36XDISTI(1|.SETUP(I).I- 1 .NG 
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AA-RHO,BETAC“2-(CG(l*l>'N(I.I)>"2-RHO,GRAV<BETAD-DY 
AB-AR0LL<I*1  )■  AA/(20*PERIOD) 


IF(AROLKI).LT.0,0)  AROLLID-OO 


[F(ABS|SBrUPNEW.SETUPOLDVSBIUPoS.GT.EPS)TOEN 

SETUPOLO-(SETUPNEW»SETUPOLDW.O 

SETUP(t)-SETUPNEW 
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CO-SQRKGRAV 


INTERPOLATE  PROFILES 


CALL  1NTEP(DIST.DPT.YA,H!C.IMAX) 
DPTKO-HIC 


IF(X'(IC).GE.XC)THEN 

DEN-X(IC)-X(ICM) 
AI-(XC-X(ICM>yDEN 
YC-A I ■Y(IQ*A2"Y(ICM) 


SUBROUTINE  RU 
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B<IF(DITH(M  ).LT.O.O)THEN 


! COEFFICENTS  FOR  CONTINUITY  EQUATION 


RETURN 


SUBROUTINE  FSWEEP(A.B.C,D.E,F^S.BS,CSJ>S.ES.FS,Q,NNN.IBC,IEC) 
IMPLIOTNONE 

INTEGER  ::  NNN,ILAST,I,!C,ICM.IBC,IEC 

DOUBLE  PRECISION.D1MENSION  (NNN)::  A,B,C.D,E,F,AS,BS.CS,DS,ES.FS,Q 
DOUBLE  PRECISION  = DEN1.DEN2 
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SUBROUTINE  FOR  RANDOM  NUMBER  GENERATION 
ADOPTED  FROM  NUMERICAL  RECIPE  OF  FORTRAN 


FUNCTION  ranO{idum) 

IMPLICIT  NONE 

INTEGER  idum,IA,lM,IQ.IR,MASK 
INTEGER  k 

DOUBLE  PRECISION  rani), AM 

PARAMETER  <lA-16807,IM=2l47483647,AM-l4M.IQ-l27773,IR-2 


APPENDIX  B 

SIMULATION  RESULTS  FOR  OCEAN  CITY,  MARYLAND 
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